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ABSTRACT 


A  detailed  field  and  laboratory  investigation  was  made  of  a  paisa 
complex  near  Fort  Simpson,  Northwest  Territories,  in  order  to  charac¬ 
terize  and  determine  the  origin  of  the  landform.  Permafrost  is  dis¬ 
continuous  in  the  area,  occurring  mainly  in  peatlands  that  have  dev¬ 
eloped  over  glacio -lacustrine  sediments.  Regionally  the  development 
of  permafrost  has  been  relatively  recent,  probably  not  occurring  before 
3500  years  Before  Present.  Paleoclimatic  information,  obtained  by  oxy¬ 
gen  isotope  analysis  of  peat  cellulose,  indicates  that  past  tempera¬ 
tures  have  not  been  an  impediment  to  permafrost  aggradation.  There  is 
no  evidence  for  temperatures  warmer  than  present  day  values  over  the 
past  10,000  years.  The  absence  of  permafrost  in  the  early  Holocene  is 
attributed  to  the  high  water  levels  in  the  peatlands.  Stratigraphic 
investigations  show  the  peatland  has  evolved  through  five  distinct 
vegetative  phases,  each  phase  being  attributable  to  increasingly  effic¬ 
ient  peatland  drainage.  It  was  not  until  water  table  divides  formed, 
creating  local  ombrotrophic  conditions,  that  extensive  Sphagnum  bogs 
developed.  As  water  levels  in  the  bogs  dropped,  possibly  aided  by  a 
period  of  intense  dessication,  the  unique  thermal  properties  of  Sphag¬ 
num  led  to  permafrost  development.  Isotopic  evidence  indicates  there 
has  been  continuous  upward  moisture  migration  through  the  permafrost 
ever  since  its  initial  aggradation. 

Stratigraphic  evidence  indicates  that  permafrost  became  areally 
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extensive,  but  not  continuous,  covering  at  least  5  km  in  the  study 
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area.  The  greatest  development  of  segregated  ice  occurred  along  the 
former  margins  of  the  glacial  lakes,  especially  in  areas  with  slightly 
elevated  mineral  soil.  With  climatic  warming,  beginning  approximately 
300  years  B.P.,  rapid  thermal  degradation  of  the  permafrost  occurred. 
Over  the  last  30  years,  the  areal  extent  of  frozen  ground  has  decreased 
by  approximately  25%>,  There  is  no  evidence  for  cyclical  aggradational- 
degradational  periods. 

The  paisa  complex  i^  not  a  frost  heave  feature,  but  rather  is  a 
relic  thermo-erosional  landform.  Melting  of  the  segregated  ice  in  the 
mineral  soil  only  accounts  for  approximately  one-half  of  the  observed 
relief  in  a  fen-palsa  peatland.  The  remaining  relief  is  due  to  normal 
consolidation  of  the  uplifted  peat  as  it  thaws.  The  end  result  of  a 
melted  paisa  is  a  shallow  circular  pond. 
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CHAPTER  1:  INTRODUCTION 
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1.1  Objectives 


Permafrost  is  defined  as  earth  materials  having  a  temperature  be¬ 
low  0°  G.  for  two  or  more  years  (Brown,  1963) .  Brown  (1968)  has  de¬ 
lineated  various  zones  of  permafrost  in  Canada,  on  the  basis  of  the 
distribution  and  frequency  of  perenially  frozen  ground.  In  the  southern 
part  of  the  Discontinuous  Zone,  permafrost  occurs  almost  exclusively 
in  peatlands. 

The  most  common  permafrost  body  in  a  subarctic  peatland  is  the 
paisa,  a  circular  to  elongate  mound  of  peat  and  mineral  soil  with  a 
perenially  frozen  core.  Individual  paisas  are  generally  less  than 
100  meters  across  and  1  to  7  meters  in  height  (Sjors,  19 6l;  Forsgren, 
1968),  but  may  coalesce,  forming  a  paisa  complex  greater  than  1  km  in 
area. 


Past  research  emphasis  on  paisas  has  been  on  their  morphology  and 
plant  ecology,  with  very  little  work  on  the  internal  structure,  ground 
ice  characteristics,  or  thermal  and  hydrological  regimes.  An  assumption 
often  made  in  previous  studies  is  that  the  gross  peatland  environment 
was  very  similar  in  the  past.  Local  edaphic  factors  were  therefore 
assumed  to  be  the  controlling  factors.  This  has  naturally  led  to  the 
notion  of  localized  and  cyclical  permafrost  dynamics  in  peatlands. 

The  hypothesis  of  this  thesis  is  that  modem  peatland  environments 
may  be  significantly  different  from  those  of  the  past  because  of  region¬ 
al  changes  in  the  climatic  and  hydrological  regimes.  To  understand 
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permafrost  dynamics,  in  subarctic  peatlands,  not  only  the  modem 
terrain  characteristics  but  also  the  history  of  the  peatland  environ¬ 
ment  must  be  established.  The  physical  and  chemical  characteristics 
of  the  peat  and  ground  ice  may  contain  information  on  the  peatland 
paleo environment . 

The  objectives  of  this  thesis  are: 

1.  To  characterize  the  terrain  of  a  paisa  complex  peatland; 

2.  To  investigate  the  genesis  of  the  paisa  complex  within  the  context 
of  the  evolution  of  the  entire  peatland; 

3.  To  determine  the  role  of  regional  climatic  fluctuations  in  perma¬ 
frost  genesis  in  the  Upper  Mackenzie  Valley; 

4.  To  determine  whether  the  formation  of  paisas  is  a  dynamic  and  cyclic 
process  or  whether  paisas  are  relic  and  erosive  landforms; 

5.  To  develop  a  conceptual  model  of  permafrost  aggradation  and 
degradation  in  a  subarctic  peatland. 

1.2  Location  of  Study  Site 

The  paisa  complex  is  located  approximately  40  km  southeast  of 
Fort  Simpson,  N.W.T.  about  midway  between  the  Mackenzie  and  Liard 
Rivers  (6l°  28*  N,  120°  55'  Ej  Fig.  l). 

The  study  site  is  located  on  a  large  speckled  bog  (Reid,  1977) 
paisa  complex  that  is  approximately  2000  m  by  250  m  in  size.  A 


. 


FIGURE  1  :  Location  of  the  study  site 


UNFROZEN, UNSATURATED 
COLLAPSED  AREA 
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representative  section  of  the  paisa  complex,  approximately  200  m  by 
100  m  was  selected  for  detailed  study.  A  sketch  map  of  the  site  is 
illustrated  in  figure  2. 


1.3  Field  Methods 

The  preservation  of  frozen  cores  is  essential  if  chemical  and 
isotopic  analysis  of  the  ground  ice  is  planned.  The  winter  drilling 
program  was  carried  oat  with  an  all  terrain  vehicle  mounted  drill, 
similar  to  that  described  by  J,  Veillette  (1975) >  capable  of  obtaining 
9  cm  diameter  core  using  a  standard  permafrost  coring  barrel.  The 
cores  were  sealed  in  plastic  film  then  packed  with  snow  in  foam  in¬ 
sulated  crates.  The  cores,  shipped  by  air  to  Edmonton,  were  exposed 
to  above  freezing  temperatures  for  less  than  8  hours  and  suffered  no 
deterioration.  The  drill  was  also  used  to  install  seven  pneumatic 
groundwater  piezometers  and  three  thermistor  cables.  Descriptions  of 
these  systems  are  included  in  chapter  two. 

Summer  field  work  took  place  over  two  months  during  1975  and  1976. 

Frozen  peat  stratigraphy  and  permafrost  depths  were  determined  by 

coring  with  a  Hoffer  probe.  Unfrozen  peat  samples  were  obtained  with 

a  standard  5  cm  diameter  piston  corer.  Permafrost  configuration  was 
extrapolated  between  cores  by  means  of  electroresistivity  profiling. 

Vegetation  sampling  was  by  a  random  point  method  developed  by  Braun - 
Blanquet  (Kershaw,  1969). 
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1.4  Previous  Theories  of  Paisa  Development 


The  role  of  vegetation,  especially  mosses  of  the  genus  Sphagnum, 
in  creating  a  micro-climate  conducive  to  permafrost  aggradation  is 
central  to  most  theories  of  paisa  formation.  In  North  America,  the 
role  of  Black  Spruce  in  establishing  a  local  negative  heat  balance 
has  been  emphasized  (Sjors,  1963;  Brown,  1968;  Zoltai,  1971;  and 
Zoltai  and  Tamocai,  1975)*  By  shading  underlying  peat  from  snow 
accumulation,  there  is  deeper  penetration  of  the  winter  frost  front. 

As  summer  melting  is  reduced  by  the  insulative  properties  of  the  peat 
and  the  shading  provided  by  the  spruce,  small  ice  lenses  may  remain 
throughout  the  year.  The  coalescing  of  these  frozen  hummocks  creates 
a  peat  plateau.  The  process  is  believed  to  initiate  with  a  single 
Black  Spruce  seedling  growing  on  a  small,  individual,  slightly  raised 
and  drier  hummock  of  Sphagnum. 

In  Fennoscandia  and  Labrador-Ungava,  paisas  are  found  north  of 
the  tree  line,  apparently  developing  without  the  aid  of  Black  Spruce. 
Under  treeless  conditions,  wind  erosion  of  the  protective  snow  cover 
on  exposed  paisa  summits  allows  deeper  propagation  of  the  frost  front 
(Salmi,  1970).  While  this  theory  explains  the  maintenance  of  the 
feature,  it  does  not  explain  the  initiation. 

The  presence  of  lichen  communities  on  paisas  has  led  to  the  hyp¬ 
othesis  that  albedo  is  an  important  factor  in  maintaining  a  negative 
heat  balance.  Lundqvist  (1951)  and  Sjors  (l96l)  have  proposed  that 
the  high  albedo  of  Sphagnum  fuscum  and  Gladonia  alpestris  is  suificient 


in  northern  regions  to  cause  development  of  an  ice  core.  Railton 
(1973)  in  a  heat  "balance  study  of  a  paisa  in  the  Hudson  Bay  Lowland 
discovered  no  evidence  for  this  and  found  variations  in  heat  balance 
from  evapo transpiration  to  be  of  more  importance. 

A  major  controversy  in  paisa  formation  is  whether  the  formation 
of  paisas  is  a  cyclic  process,  as  first  advocated  by  Lundqvist  (1953.) 
or  whether  paisas  are  relic  features,  formed  in  a  past  period  of  cli¬ 
matic  deterioration.  Lundqvist' s  theory  of  cyclical  aggradation  and 
degradation  has  been  advanced  in  North  America  by  Sjors  (1963),  Brown 
(1968),  Zoltai  (1971),  Zoltai  and  Tarnocai  (1971  and  1975) >  Payette 
and  Samson  (1976),  Railton  (1973) »  and  Reid  (1977).  Zoltai  and  Tarn¬ 
ocai  (1975)  have  applied  the  terms  incipient,  young,  mature  and  over¬ 
mature  to  the  developmental  stages  of  paisas  and  peat  plateaus.  The 
incipient  and  young  peat  plateaus  are  considered  to  be  an  aggrading 
stage.  The  mature  plateau  is  in  a  state  of  equilibrium  which  upon 
reaching  the  degrading  stage  is  referred  to  as  overmature.  In  con¬ 
trast  to  this  theory,  Vorren  (1972)  and  Sollid  and  Sorbel  (197^) 
support  the  idea  that  paisas  are  relic  landscape  features  which  are 
out  of  phase  with  present  climatic  conditions,  and  that  the  ground 
ice  is  protected  only  by  the  insulating  peat  layer. 

Further  controversy  in  paisa  development  is  that  of  a  dynamic 
versus  an  erosive  origin.  The  dynamic  theory,  first  postulated  by 
Lundqvist  (l95l),  involves  the  progressive  growth  of  a  small  Sphagnum 
frost  hummock  to  a  mature  paisa.  The  previously  mentioned  process  in¬ 
volving  a  Black  Spruce  initiated  micro-climate,  is  a  dynamic  process. 
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This  theory  has  received  wide  spread  support  both  in  Europe  and  in 
North  America  (Sjors,  19 6l;  Zoltai,  1971;  Zoltai  and  Tarnocai,  1975)* 

A  theory  first  suggested  by  Cajuander  (1913)  is  that  paisas  might 
be  erosional  features  in  the  process  of  degrading.  The  region  where 
present  day  paisas  are  found  is  thought  to  have  once  been  covered  by 
a  broad  area  of  polygonal  frozen  ground.  Climatic  warming  has  pro¬ 
duced  ground  ice  degradation  and  thermokarsting.  Paisas  are  interp¬ 
reted  as  erosional  remnants  of  the  once  more  extensive  permafrost. 

More  recently,  Vorren  (1972)  has  advocated  an  erosive  origin  for 
paisas  in  the  northern  part  of  Norway. 
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2.1  Climate 


The  study  site  is  within  the  Sub-Boreal  Forest  Climatic  Zone 
(Burns,  1973)*  Mean  daily  temperature  at  Fort  Simpson,  40  km  to  the 
north-west,  is  -4.1°  C.j  the  mean  July  temperature  is  16°  C.;  and  a 
mean  temperature  of  -29°  C.  is  recorded  for  January.  The  number  of 
freezing  degree-days  (below  32°  F.)  is  6300  and  the  number  of  thawing 
degree-days  (above  32°  F.)  is  34.50,  Precipitation  averages  34.5  cm 
of  which  40%  falls  as  snow. 

2.2  Glacial  History  and  Geology 

The  area  was  invaded  at  least  twice  from  the  east  by  Laurentide 
Ice  (Rutter  et  al,  1974).  The  first  advance  is  recorded  by  a  grey- 
black  basal  stony  till  with  erratics  found  up  to  an  elevation  of 
1500  m  in  the  Mackenzie  Mountains  to  the  west.  The  second,  lesser, 
advance  believed  to  be  Classical  Wisconsin,  deposited  a  light  grey- 
brown  stony  till  (Rutter  et  al,  1973)* 

Retreat  of  the  ice  sheet  was  northeastward,  damming  natural 
drainage  outlets  and  creating  large  glacial  lakes  along  the  Mackenzie 
and  Liard  River  Valleys.  Glacial  Lake  McConnell  was  the  largest  of 
these  lakes,  remaining  until  at  least  11,500  years  B.P.  (Rutter  et  al, 
1973) •  Large  fluctuations  in  lake  levels  created  well  defined  beach 
ridges . 


' 


The  study  site  is  located  approximately  250  meters  north  of  one 
of  the  prominent  beach  ridges.  The  beaches,  composed  of  fine  to  medium 
grained  sands  with  occasional  gravel,  are  up  to  6  m  in  height.  To  the 
north,  is  a  broad  silty  clay  glaciolacustrine  plain  capped  uniformly 
with  organic  deposits.  Occasional  sand  dunes  provide  the  only  relief. 
Surficial  sediment  is  underlain  by  Lower  Cretaceous  shales  of  the  Ft. 
St.  John  group  and  by  Upper  Devonian  cryptocrystalline  and  clastic 
limestones  of  the  Dunvegan  formation. 


2.3  Landscape  Position 

The  areal  distribution  of  six  peatland  types  is  superimposed  on 
the  surficial  geology  map  of  the  area  (Rutter,  et  al,  1973?  Fig.  3)* 

The  peatland  types  are: 

1.  fen 

2.  fen  with  scattered  peat  plateaus 

3 .  peat  plateau 

4.  'speckled  bog' 

5.  'speckled  bog'  paisa  complex 

6.  collapse  scars  and  ponds. 

They  are  subdivided  on  the  basis  of  ground  ice  development,  but  mapped 
by  air  photo  identification  of  associated  surface  characteristics  (see 
legend,  Figure  3)» 


There  are  definite  physiographic-peatland  type  relationships  ( see 
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Figure  3  :  Distribution  of  peatland  types  and  associated  beach  deposits 
in  the  Ft.  Simpson  area. 
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Figure  3) •  Of  particular  interest  is  the  association  of  the  paisa  com¬ 
plexes,  representing  the  greatest  segregational  ice  development,  with 
the  inner  edges  of  former  glacial  lake  margins.  The  two  areas  of  best 
paisa  development  are  the  study  area  at  the  southern  edge  of  the  lake 
basin  and  near  Antoine  Lake  at  the  western  edge. 

Collapse  scars  and  non-oriented  lakes  (peatland  type  C,  Figure  3) 
are  intimately  associated  with,  often  completely  surrounding,  the  paisa 
complexes.  The  central  areas  of  the  lake  basin  contains  lakes,  un¬ 
frozen  fens  and  scattered  peat  plateaus.  The  peripheral  edges  of  the 
lake  margins  are  covered  in  shallow  'peat  deposits  and  modem  day  sedge 
fens.  The  permafrost  in  these  areas  is  scattered  and  very  shallow. 

The  till  uplands  are  either  permafrost  free  in  the  well  drained 
sites  or  contain  large  ’speckled  bog’  peat  plateaus  in  low  lying  de- 
pressional  areas. 

It  must  be  emphasized  that  this  is  a  very  generalized  zonation. 
Individual  paisas  cannot  be  identified  at  this  scale  o*  mapping;  there* 
fore  they  may  occur  under  certain  topographical  and  hydrological  con¬ 
ditions  in  the  more  central  parts  of  the  basin.  Individual  paisas  have 
been  observed  on  the  till  plains,  especially  where  there  is  abundant 
water.  Paisa  complexes  greater  than  100  m  in  diameter  are,  however, 
restricted  to  the  margins  of  the  former  glacial  lake. 


■ 


2.4  Morphology 


The  study  site  was  selected  to  represent  all  major  geomorphic 
features  of  the  permafrost  body.  Morphology,  vegetation  and  ground 
ice  descriptions  are  applicable  to  the  entire  paisa  complex.  The  open 
fen  Surrounding  the  paisa  is  distinctly  different  however,  and  is  de¬ 
scribed  separately. 

The  paisa  complex  is  approximately  2000  m  long.  Its  width  at 
the  broadest  point  is  about  300  m  "but  narrows  to  100  m  at  its  eastern 
end.  The  total  area  is  about  4  hectares. 

Sharp  collapsing  edges  with  an  average  height  of  2.5  meters,  mark 
the  edge  of  the  paisa.  The  edges  are  discontinuous  however;  with  more 
gently  sloping  peat  banks  common. 

Three  transects  (Figure  2)  were  made  to  determine  the  surface 
morphology,  the  active  layer  depths,  the  depths  of  peat  and  total 
ground  ice  thickness.  The  transects  are  illustrated  in  three  cross- 
sections  (Figures  4,  5  and  6). 

Surface  elevations  range  from  2.5  to  4.7  meters  above  the  sur¬ 
rounding  fen,  averaging  3*0  meters.  There  is  considerable  local  re¬ 
lief  however,  created  mainly  by  interior  thawed  Sphagnum  bogs  and  by 
interior  ' thermokarst '  sinkholes.  The  sinkholes  have  near  vertical 
banks  and  show  signs  of  active  collapse.  The  Sphagnum  bogs  have  more 
gently  sloping  banks  and  are  usually  less  than  1  meter  below  the  sur- 
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FIGURE  *3  «  TOPOGRAPHIC  PROFILE,  GROUND  ICE  CONFIGURATION,  GEOPHYSICAL  RESISTIVITY  RESPONSE. 

PEAT  STRATIGRAPHY  AND  VEGETATION  DISTRIBUTION  FOR  TRANSECT  S — S' 
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FIGURE  e  .  TOPOGRAPHIC  PROFILE  .  GROUND  ICE  CONFIGURATION  .  GEOPHYSICAL  RESISTIVITY 

and  vegetation  distribution  for  transect  n  n  . 


22 


face  of  the  frozen  peat. 

Paisa  heights  appear,  from  the  cross-sections,  to  correlate  with 
the  total  thickness  of  permafrost.  Where  permafrost  is  thicker,  ground 
elevations  are  higher. 

Distinct  marginal  ridges  are  elevated  approximately  0.5  to  1.0 
meters  above  the  interior  regions  of  the  paisa  complex.  The  ridges 
are  discontinuous,  usually  associated  with  steep  collapsing  edges. 

Section  E-E'  (Figure  4)  contains  an  unfrozen  unsaturated  collapsed 
bog.  The  peat  stratigraphy  indicates  a  former  saturated  Sphagnum  bog 
that  has  drained  due  to  melting  of  the  adjacent,  damming  permafrost. 

The  surface  elevation  is  locally  depressed  to  nearly  that  of  the  sur¬ 
rounding  fen.  The  feature  is  marked  by  sharp  collapsed  edges  on  3 
sides  and  a  'hinge-like'  connection  on  the  fourth  side. 

Active  layer  depths  are  strongly  controlled  by  vegetation.  Indiv¬ 
idual  paisa  summits,  which  are  sparsely  treed,  have  the  greatest  active 
layers,  averaging  85  cm  in  late  August.  Active  layers  average  55-6 0  cm 
elsewhere,  except  in  areas  of  dense  Black  Spruce-  feathermoss  com¬ 
munities  where  the  mean  depth  is  only  35  cm. 

The  surface  of  the  paisa  complex  is  extensively  dissected  by  cracks 
extending  the  full  thickness  of  the  active  layer.  Crack  densities  and 
dimensions  are  greatest  on  paisa  summits,  especially  near  collapsing 
edges  where  a  strong  orientation  parallel  to  the  edge  is  observed.  The 
moister  regions  of  the  paisa  complex,  especially  in  the  dense  Black 
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Spruce-f eathermoss  communities,  are  less  dissected  tending  to  have 
cracks  less  than  7  cm  in  width  compared  to  widths  of  over  70  cm  on 
paisa  summits.  The  role  of  surface  cracking  on  ground  ice  degra¬ 
dation  is  discussed  in  chapter  four. 

The  unfrozen  peatland  surrounding  the.  paisa  complex  is  an  open 
sedge  fen.  The  fen  consists  of  small  randomly  oriented  ponds,  often 
greater  than  3  m  in  depth,  separated  hy  floating  mats  of  sedges.  The 
fen  area  on  the  beach  ridge  side  (south  side)  of  the  paisa  complex 
appears  to  have  recently  been  much  drier,  probably  a  flat  fen  or 
horizontal  fen  (Tamocai,  1970  classification).  The  area  once  sup¬ 
ported  a  birch  grove,  however  renewed  paludification,  probably  due  to 
beaver  dams,  has  reestablished  the  open  fen.  Aerial  photography  shows 
the  area  has  been  flooded  since  at  least  1946.  The  thickness  of  peat 
in  the  fen  is  less  than  on  the  paisa  complex,  averaging  only  1.5 
meters,  compared  to  the  3«0  meter  average  on  the  paisas. 

Peat  plateaus  are  scattered  throughout  the  fen.  There  is  a  great 
variation  in  size,  from  a  maximum  of  1  hectare  to  a  minimum  of  only  a 
few  tens  of  square  meters,  but  averaging  approximately  400  square 
meters.  The  surface  elevation  of  the  peat  plateaus  is  much  less  than 
the  paisa  complex,  averaging  1  m  above  the  surrounding  fen.  Bank  edges 
lack  a  steep  collapsing  bank.  Thermokarsting  is  evident  however  from 
the  numerous  dead  trees  leaning  toward  the  fen. 


# 


% 


2 .3  Vegetation 


Vegetation  on  the  paisa  complex  and  on  surrounding  fen  was  sampled 
using  a  random  point  sample  technique.  Percent  cover  of  each  species 
was  visually  estimated,  then  recorded  on  a  five  point  scale  modified 
from  Braun -Blanquet  (Kershaw,  1969).  The  main  peatland  types  with 
their  respective  vegetative  communities  are  detailed  in  Table  1. 


2.6  Stratigraphy 

The  peat  stratigraphy  of  the  paisa  complex  was  determined  by 

macroscopic  analysis  of  46  cores.  Depth  profiles  along  two  transects 

are  illustrated  in  Figures  4  and  3*  The  classification  of  the  peat 

materials  is  similar  to  that  used  by  Zoltai  and  Tarnocai  (1975)*  The 

peat  types  and  their  modem  growth  environments  are  summarized  in 
Table  2. 

The  peat  strata  record  a  sequence  of  six  broad  vegetative  com¬ 
munities  in  the  following  order  from  the  bottom  of  the  core  to  the  top: 

1.  detrital  organic  in  a  silt  matrix 

2.  aquatic 

3.  brown  moss-sedge  fens 

4.  sedge-shrub  communities 
3 .  Sphagnum  bogs 

6.  Black  Spruce-feathermoss-Cladonia  communities. 

Each  of  these  communities  implies  a  different  wetland  environment. 


PORTLAND  TYPES  AND  VEGETATION  COMMUNITIES 
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Modem  correlatives  of  the  vegetation  sequence  implies  an  environmental 
transition  through: 

1.  large  deep  lake 

2.  numerous  shallow  lakes  and  ponds 

3.  open  sedge  fen 

4.  sedge  and  shrub  fen  or  paludificated  shrub  forest  or  fen  ridge 
complex 

5*  ombro trophic  Sphagnum  bogs 
6.  raised  frozen  peat  plateau. 

Variation  is  found  between  cores  however,  implying  that  local 
edaphic  factors  have  played  a  part  in  determining  the  local  stratig¬ 
raphy.  For  example,  in  some  cores,  the  woody  fen  layer  is  absent  or 
is  succeeded  by  a  brown  moss-sedge  peat.  The  above  scheme  does,  how¬ 
ever,  portray  the  average  environmental  sequence. 

The  mineral  soil  stratigraphy  reflects  the  proximal  position  of 
the  paisas  in  the  former  glacio lacustrine  basin.  The  contact  between 
peat  and  mineral  soil  is  sharp.  Typically,  the  underlying  sediment 
is  a  silty  sand  for  the  uppermost  .3  to  1.0  m,  grading  to  a  uniform, 
massive,  blue-grey  silty  clay.  The  grain  size  distribution  of  a 
composite  sample  of  the  clay  is  64%  clay,  32%  silt  and  4%  sand. 
Occasional  lenses  of  fine  to  medium  grained  sand  randomly  intersect 
the  silty  clay.  Large  water  yields  were  encountered  when  unfrozen 
sand  lenses  were  drilled,  indicating  the  lenses  are  fairly  continuous. 
Glacial  till  was  not  encountered  up  to  a  depth  of  11  meters.  Nearby 
shot  hole  data  indicate  approximately  15-20  m  of  glaciolacustrine  clays 
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overlie  approximately  3  to  10  m  of  till. 


2.7  Permafrost  Configuration 

The  configuration  of  the  permafrost  body  is  illustrated  in  three 
cross-sections  (Figures  4,  5  and  6).  The  permafrost  boundaries  were 
determined  through  interpretation  of  an  electromagnetic  resistivity 
survey  completed  with  a  Geonics  EM  31  instrument.  The  predicted 
depth  of  permafrost  correlates  well  with  the  depth  measured  in  44 
boreholes.  Boundaries  between  frozen  and  unfrozen  zones  are  sharp 
and  well  defined. 


2.8  Ground  Thermal  Regime 

Permanent  thermistor  strings  were  installed  to  measure  the 
thermal  structure  of  the  paisa  and  the  adjacent  fen  and  collapse  bog. 
Random  ground  temperatures  were  also  measured  with  a  needle  probe 
thermistor.  The  method  of  thermistor  construction  is  detailed  in 
Appendix  1. 

The  vertical  temperature  profiles  of  the  frozen  paisa,  to  a 
depth  of  15  m,  and  the  adjacent  unfrozen  bog,  to  a  depth  of  9 
are  illustrated  in  Figures  7  and  8  respectively.  Temperatures  were 
recorded  in  mid-March  and  late  August,  hence  should  reasonably  rep¬ 
resent  the  annual  soil  temperature  variation  (see  Brown,  1966) . 
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FIGURE  B  ■  TEMPERATURE  PROFILE  IN  UNFROZEN 

PEAT  BOG  WITHIN  PALSA  COMPLEX. 
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Thermal  characteristics  of  the  paisa  and  unfrozen  hog  are: 


Frozen  Paisa 

Unfrozen  Bog 

1. 

Mean  annual  air  temperature 
(obtained  from  Ft.  Simpson) 

-  4.1°  G 

-  4.1°  G 

2. 

Mean  annual  surface  temperature 

-  2.1°  G 

+  .9°  C 

3. 

Temperature  gradient 

+  .18°  C/m 

+  .05°  C/m 

4. 

Depth  of  annual  temperature 

variation 

3.4  m 

1.6  m 

The  mean  annual  surface  temperature  of  the  paisa  is  significantly 
colder  than  mineral  soil  surface  temperatures  at  Ft.  Simpson,  of  +  0.8 
to  +  1.6°  C.  (Bums,  1974).  This  is  consistent  with  other  studies  in 
peatlands  (Williams,  1968). 

The  temperature  gradient  in  the  paisa  is  unusually  steep  compared 
with  other  measurements  made  in  the  Ft.  Simpson  area  (Judge,  1974). 

The  gradient  in  shallow  permafrost  is  very  sensitive  to  changes  in 
mean  winter  temperatures  because  of  the  high  thermal  conductivity  and 
Sphagnum  cover.  Sphagnum  promotes  rapid  outward  heat  flow  during 
winter  because  its  thermal  conductivity  when  frozen  is  ten  times  that 
when  dry  and  unfrozen  (JumiJkis,  1977).  The  recent  decline  in  mean 
winter  temperatures  in  the  Ft.  Simpson  area  over  the  past  eight  years 
(Bums,  1974;  Atmospheric  Environment  Service  Records  1974-  1977)  has 
probably  created  the  steep  temperature  gradient  observed  in  the  paisa. 
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2.9  Ground  Ice  Characteristics 

i 

2.9I  Description 

The  frozen  peat  and  mineral  soil  contain  a  variety  of  ground 
ice  types.  Fourteen  cores,  each  9  cm  in  diameter,  were  examined  in 
the  field,  to  determine  ground  ice  variability  across  the  paisa.  The 
permafrost  ms  classified  as  pore  ice,  vein  ice,  reticulate  vein  ice, 
or  as  ice  lenses  of  a  specified  thickness  (Figure  9  ).  Lenses  thicker 
than  8  cm  are  illustrated  individually. 

Pore  ice  is  primarily  restricted  to  the  active  layer.  The  active 

layer  is  usually  within  the  unsaturated  fibric  Sphagnum  horizon,  with 

little  moisture  available  for  ice  lens  development.  Crystal  sizes, 

2 

averaging  1  mm  or  less,  were  equidimensional  throughout  the  entire- 

active  layer.  In  one  core,  slightly  larger  crystals  with  small 

spherical  bubbles  had  developed  within  a  forest  peat  layer.  The  small 

crystals  and  small  spherical  bubbles  represent  rapid  freezing  in  an 
undersaturated  medium  (Gell,  1976). 

Most  of  the  peat  column  contains  numerous,  closely  spaced,  friable 
lenses  less  than  2  mm  in  thickness.  The  lenses  are  oriented  parallel 
to  the  inferred  heat  flow  direction.  Bubbles  are  numerous  and  appear 
spherical.  Water  content  averages  800%  by  weight  and  46%  by  volume 
which  is  approximately  the  same  as  the  adjacent  unfrozen  Sphagnum  bogs. 

The  basal  0.5  m  of  peat  often  contains  numerous  ice  lenses  up  to 
10  cm  in  thickness.  Crystals  are  only  slightly  elongate,  averaging  15 
to  20  mm  in  length  and  15  mm  in  width.  The  basal  peat  is  considerably 
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more  decomposed  (H  =  6,  von  Post  Scale)  than  the  upper  portions  of  the 
peat  column.  The  greater  amounts  of  colloidal  sized  particles  and 
slowly  moving  freezing  fronts  have  combined  to  encourage  segregational 
ice  lensing.  Ice  lensing  results  from  the  nucleation,  growth  and 
coalescence  of  ice  crystals  in  layers  normal  to  the  direction  of  heat 
extraction  in  freezing  soil  (Gell,  1976).  .Subtle  variations  of  temp¬ 
erature  at  the  freezing  plane  are  believed  to  determine  whether  seg¬ 
regational  ice  grows  in  the  field  situation  (Mackay,  197*0  •  This  may 
explain  why  lensing  in  the  lowermost  peat  is  highly  variable,  creating 
lenses  up  to  10  cm  in  some  cores  but  not  creating  any  lenses  in  other 
cores. 

A  characteristic  of  most  cores  was  a  large  ice  lens  at  the  peat 
mineral  soil  interface.  The  lens  was  up  to  85  cm  in  thickness,  typ¬ 
ically  contained  peat  inclusions  at  the  top  but  was  clear  for  the  bottom 
75  Per  cent.  Crystals  were  large,  some  up  to  8  cm  long  and  2.5  cm  wide, 
witn  the  c-axis  strongly  oriented  in  the  vertical  direction.  Although 
few  were  seen,  bubbles  were  slightly  elongate  parallel  to  the  c-axis. 
These  characteristics  are  typical  of  slow  crystal  growth  (Anderson  and 
Morgenstern,  1973)* 

The  mineral  soil  is  characterized  by  numerous  discrete  ice  lenses 
separated  by  overconsolidated  mineral  soil  containing  reticulate  ice 
veins.  In  some  cores,  such  as  N9  and  S56,  the  ice  lenses  are  spaced 
nearly  evenly  throughout  the  core.  In  other  cores  however,  such  as 
S78,  S100  and  E3»  lenses  are  concentrated  in  the  upper  portion  of  the 
core  with  reticulate  vein  ice  dominating  the  lower  half  of  the  frozen 


. 
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mineral  soil.  This  is  the  exact  opposite  of  what  would  he  expected 
with  a  continually  slowing  freezing  front  with  depth  (Mackay,  197*0  • 
Similar  profiles  have  been  observed  by  Gell  (1976),  in  pingo  ice,  who 
ascribed  them  to  an  increase  of  freezing  rates  due  to  pingo  uplift 
and  exposure  of  the  surface  to  cold  air  temperatures .  A  similar 
mechanism  would  account  for  the  basal  reticulate  vein  ice  in  the  paisa 
cores,  as  could  a  series  of  very  cold  winters  with  little  snowfall. 

Many  of  the  cores  (S78,  S56,  SlOOa,  Ell  and  N9)  have  a  large  ice 
lens  averaging  .4m  thick,  but  up  to  .75  m  thick,  just  above  the  base  of 
the  permafrost.  This  lens,  typically  thicker  than  the  other  mineral 
soil  lenses,  is  often  bounded  by  long  sections  of  reticulate  vein  ice. 

Crystal  structures  in  the  segregational  vein  ice  are  well  developed 
with  an  assortment  of  sizes.  Maximum  observed  crystal  dimensions  were 
70  mm  in  length  and  15  mm  in  width;  c-axis  orientations  were  all 
parallel  to  the  core.  Macroscopic  bubbles  were  rare  and  slightly  elon¬ 
gate  where  observed. 


2.92  Ionic  Chemistry 

The  major  ion  chemistry  of  the  paisa  ground  ice  can  aid  in 
determining  the  hydrological  history  of  the  site.  Vertical  variability 
in  conductivity  of  permafrost  cores  may  represent  periods  of  thaw  and 
surface  leaching  (O'Sullivan,  1963;  Pewe  and  Sellman,  1972);  while 
horizontal  variations  may  reflect  differences  in  water  sources  (Brown, 


■ 


1963;  Allen,  1976). 


The  hydrochemistry  of  eight  potential  water  sources  for  the  ground 
ice  was  determined.  In  addition,  six  permafrost  cores  were  analyzed 
for  major  ions,  pH  and  conductivity.  The  major  cations:  calcium, 
magnesium,  sodium  and  potassium  were  analyzed  by  atomic  absorption 
spectroscopy.  Chloride  was  determined  by  mercuric  nitrate  titration; 
and  bicarbonate  by  potentiometric  titration  with  sulphuric  acid.  Sul¬ 
phate  ion  concentration  was  determined  gravimetrically  using  a  barium 
chlorate  reagent  for  the  mineral  soil  ice  and  by  difference  for  the 
peat  ice.  As  a  check  on  the  validity  of  the  analysis,  cation  and  an¬ 
ion  balances  were  calculated  for  each  sample.  In  general,  the  balance, 

expressed  as  cations,  ^  anions  was  less  than  2%  for  80 %  of  the 

cations  +  anions 

samples  and  less  than  for  the  remainder.  Analytical  results  are 
tabulated  in  Appendix  Two. 

The  potential  water  sources  for  the  ground  ice  have  quite  distinct 
chemistry  (Table  3)«  Sphagnum  bog  water  is  very  similar  to  that  of 
precipitation;  conductivity  averages  40  mhos/cm  and  pH  is  approx¬ 
imately  5.0.  Water  from  the  fen  has  a  conductivity  of  250  mhos/cm, 
with  a  pH  of  7.9;  calcium  and  bicarbonate  are  the  major  ions.  Deeper 
mineral  soil  groundwater  has  moderately  high  conductivities  (350-550 
mhos/cm,  and  a  pH  of  7.4.  Calcium,  magnesium,  and  sodium  are  the 
major  cations;  bicarbonate  is  the  major  anion. 


Hydrochemistry  of  the  permafrost  cores  reflects  the  various  source 
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waters  in  some  instances,  but  is  significantly  enriched  in  other  parts 
of  the  core.  The  upper  frozen  peat  portions  of  the  cores  are  very  low 
in  dissolved  solids  (Figures  10  to  13),  similar  to  the  modem  Sphagnum 
bogs.  This  indicates  that  at  the  time  of  freezing,  the  bogs  were 
probably  ombrotrophic  as  their  surfaces  were  raised  above  the  regional 
water  table,  therefore  receiving  no  mineralized  groundwater. 

The  chemistry  of  the  lower  frozen  peat  cores  is  again  similar  to 
modem  fen  water,  except  for  a  marked  increase  in  the  S0^  ion.  As 
modem  peatland  water  and  subpermafrost  groundwater  shows  only  trace 
amounts  of  S0^  ,  the  origin  of  the  S0^  ion  is  problematic.  One  possib¬ 
ility  is  the  mineralization  of  sulphate  which  occurs  during  decompos¬ 
ition  of  organic  matter  (Barrow,  i960  and  1961).  In  non-decomposed 
plant  material,  sulphur  is  present  as  organically  bound  sulphate. 

With  decomposition,  available  carbon  is  progressively  utilized  by  soil 
microbes  until  the  carbon/sulphur  ratio  is  about  50,  after  which  sulphur 
mineralization  occurs.  Barrow  (1961)  also  found  a  large  increase  in 
inorganic  sulphur  with  drying  of  the  organic  matter.  Freezing  has  been 
shown  to  be  a  similar  physicochemical  process  to  drying;  while  ex¬ 
perimental  evidence  is  lacking,  it  would  seem  reasonable  to  predict  an 
inorganic  SO^  increase  during  freeze-thaw  cycles.  Tamocai  (197*0  bas 
found  such  an  increase  in  Ca‘+  content  of  frozen  organic  soils  which  he 
attributes  to  the  freezing  process.  The  origin  of  the  sulphate  ion  in 
the  decomposed  peat  is  further  supported  by  comparison  of  the  ground- 
water  chemistry  beneath  a  frozen  paisa  with  that  beneath  a  recently 
collapsed  paisa:  beneath  the  frozen  paisa,  sulphates  average  0*5  epm, 
while  beneath  the  collapsed  paisa,  sulphates  average  1.4  epm. 
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FIGURE  1  1  *  MAJOR  ION  CHEMISTRY  PROFILE  FOR  CORE 
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FIGURE  13  ■  MAJOR  ION  CHEMISTRY  PROFILE  FOR  CORE  S-78. 
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The  most  striking  feature  of  four  of  the  cores  is  the  strong  peak 
in  Ga  and  S0^_  displayed  at  the  peat -mineral  soil  boundary.  Part  of 
this  peak  is  probably  due  to  the  above  freeze-thaw  mechanism,  however 
the  extreme  concentration  is  most  likely  the  result  of  solute  precip¬ 
itation  as  gypsum  (CaSO^).  This  precipitation  would  only  occur  during 
a  period  of  intense  dessication.  Gypsum  crystals  have  been  so  observed 
by  F.W.  Schwartz  (personal  communication)  while  drilling  in  dry  lake 
beds  in  central  Saskatchewan.  The  timing  of  the  dessication  period, 
based  on  evidence  from  central  Alberta  (May  and  Thomson,  1977 )>  was 
probably  immediately  post-glacial  preceeding  the  build-up  of  peat. 


2.93  Environmental  Isotopes 

Both  stable  and  natural  radiogenic  isotopes  are  well  estab¬ 
lished  as  hydrologic  research  tools.  Recently,  stable  isotopes  have 
also  provided  information  on  the  thermal  history  (Michels  and  Fritz, 
1978)  and  freezing  mechanisms  (Mackay  and  Lavkulich,  1974)  of  perma¬ 
frost.  In  spite  of  only  recent  application,  the  potential  use  of  iso¬ 
topes  in  determining  past  temperatures  was  established  over  20  years 
ago  by  Weston  (1955)  and  others.  Water  was  found  to  undergo  a  maximum 
2  %  shift  in  its  01  /  01  ratio  during  freezing,  thereby  maintaining 
a  record  of  its  original  isotopic  composition,  which  in  turn  is  a 
function  of  the  air  temperature  at  the  time  of  precipitation  (see 
O’Neil,  1968  ). 


Natural  radiogenic  isotopes  have  also  seen  very  limited  use  in 
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permafrost  research,  "being  used  in  a  tritium  study  "by  Forsgren  (1964) 
in  Northern  Sweden.  Because  of  the  relatively  short  half-life  of 
tritium  (12.26  years)  its  use  as  a  dating  isotope  is  restricted.  It 
can,  however,  be  employed  as  an  indicator  of  recent  or  old  permafrost. 
This  is  because  of  the  enormous  increase  in  the  natural  tritium  content 
of  meteoric  water  that  has  occurred  since  the  detonation  of  nuclear 
devices  in  the  early  1950’s.  For  example,  precipitation  at  Ft.  Smith, 
N.W.T.  reached  a  post-bomb  maximum  of  3616  TU  in  1963,  which  has  since 
decreased  to  about  200  TU.  According  to  Simard . (1977) >  groundwater  with 
concentrations  greater  than  20  TU  has  been  contaminated  by  thermonuclear 
tritium  hence  are  of  recent  origin;  water  with  concentrations  less  than 
3  TU  are  said  to  be  free  tritium  water  introduced  into  the  flow  system 
prior  to  the  bomb  testing  period;  tritium  concentrations  between  3  TU 
and  20  TU  indicate  mixing  of  post  bomb  and  pre-bomb  precipitation. 

Oxygen,  deuterium,  and  tritium  isotopic  analyses  were  made  on  a 
permafrost  core  from  the  paisas,  as  well  as  on  several  groundwater 
samples  collected  through  piezometer.  Condensation  of  tritium-rich 
vapor  onto  the  ice  cores  was  minimized  by  coring  at  temperatures  less 
than  the  temperature  of  the  permafrost,  then  quickly  wrapping  the  frozen 
core  in  air-tight  polyvinyl  and  heavy  plastic.  Such  precautions  were 
not  necessary  for  the  stable  isotopes.  Samples  were  taken  from  both  ice 
lenses  and  fine  reticulate  vein  ice.  This  necessitated  separation  of 
soil  and  water  by  high  speed  centrifugation;  fractionation  effects 
due  to  centrifugation  are  unknown,  but  are  most  probably  negligible. 
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2.931  Tritium 

All  tritium  measurements  were  made  at  the  Weizmann  Inst¬ 
itute,  using  internal  gas  counting  preceeded  "by  electrolytic  enrich¬ 
ment.  The  analytical  results  are  listed  in  Table  Four.  Tritium  val¬ 
ues  for  the  sut-permafrost  groundwater  and  the  ice  lenses  in  the  min¬ 
eral  soil  are  essentially  equal.  The  groundwater  value  of  7.9  TU  in¬ 
dicates  that  the  water  is  a  mixture  of  pre-1950  and  post-1930  water. 

This  supports  the  theory  that  the  mineral  soil  groundwater  receives 
minimal  recharge  from  fen  surface  waters  ( see  previous  section) ,  The 
high  tritium  values  for  the  uppermost  frozen  peat  sample  indicates  in¬ 
filtration  hy  modern  precipitation.  The  depth  of  the  infiltration  ap¬ 
pears  to  he  limited  to  the  upper  0.5  m  of  frozen  peat  on  the  basis  of 
oxygen  -18  analysis. 

The  magnitude  of  the  tritium  value  of  the  ground  ice  in  the  min¬ 
eral  soil  as  well  as  the  fact  that  it  is  nearly  the  same  as  the  value 
for  the  sub-permafrost  groundwater  presents  a  paradox.  Carbon-14  dating 
of  the  paisa  indicates  an  approximate  age  of  1500  years  (see  Chapter  3) • 
Tritium  values  for  uncontaminated  ground  ice  of  this  age  should  be  less 
than  2  TU.  The  magnitude  of  the  ground  ice  value  exceeds  the  likely 
amount  of  contamination,  but  this  possibility  certainly  cannot  be  ex¬ 
cluded.  Recent  ground  ice  melting  followed  by  rapid  reheaving  so  that 
no  additional  peat  is  added  to  the  paisa,  may  also  cause  anomolous 
tritium  values.  Melting  is  not  supported  by  peat  stratigraphic  or  thaw 
consolidation  analysis,  however  (see  Chapter  5).  A  third  explanation, 
and  one  that  also  accounts  for  the  similarity  of  tritium  values  in  the 


TABLE  4 


ENVIRONMENTAL  TRITIUM  ANALYSIS 


Sample  Number  Depth  (meters)  Tritium  Value_ Location 


FC-1 

1.20 

23  -  0.3  T.U. 

Peat 

FC-2 

3.5 

7.9  -  0.3  T.U. 

Peat -mineral  soil 
interface 

FG-3 

6.5 

6.4  -  0.4  T.U. 

Mineral  soil  ice 
lenses 

FC-4 

10.0 

6.4  i  0.4  T.U. 

Sub-permaf ro  st 
groundwater 

' 
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ice  and  groundwater,  is  for  some  water  exchange  to  have  occurred  "between 
the  frozen  core  and  the  groundwater.  This  explanation  is  supported  by 
the  stable  isotope  data  which  is  elaborated  upon  in  the  following 
section. 


2.932  Stable  Isotopes 

Oxygen  isotope  ratios  were  determined  by  the  carbon  diox¬ 
ide  equilibration  technique  of  Epstein  and  Mayeda  (1953).  Deuterium 
measurements  were  made  by  reacting  the  samples  with  uranium  metal  at 
600  degrees  Celsius  to  produce  hydrogen  gas  (Eigeleisen,  1952).  The 
analytical  results  are  listed  in  Table  5>  as  well  as  plotted  as  a  depth 
profile  in  Figure  14. 

Prior  to  analyzing  the  stable  isotope  profiles,  it  is  necessary 
to  determine  whether  any  isotopic  enrichment  had  occurred  prior  to 
freezing.  Such  enrichment  can  be  caused  by  partial  evaporation  of  the 
water.  The  previous  occurrence  of  evaporation  can  be  detected  by 
measuring  both  oxygen  and  deuterium  isotope  values  for  the  same  sam¬ 
ples,  then  comparing  a  plot  of  these  values  with  the  Meteoric  Water 

18 

Line  (Craig,  1961).  A  plot  of  6  D  vs.  cf  0  values  from  the  sam- 

18 

pies  falls  on  the  regression  line:  6  D  =  8  cT~0+5  (Figure  15). 

This  equation  is  the  same  as  that  found  by  Hage  et  al  (1975)  for  weigh¬ 
ted  mean  precipitation  from  high  latitude  continental  areas.  The  slight 
deviation  from  the  meteoric  water  line  is  believed  due  to  near  equili¬ 
brium  evaporation  of  North  Pacific  sea  water  during  the  summer  (Hage 
et  al,  1975).  The  ground  ice  samples  are  therefore  representative  of 


. 


TABLE  5 


49 


£  £ 

o  o 


CD  CD 


3 

£ 

Mi 

£ 

ft 

ft 

ft 

ft 

0 

3-"-m 

3^ 

3 

3 

ccJ 

3 

•H 

CO  >5 

CO  >5 

0 

0 

0 

0 

-P 

rH 

rH 

ft 

ft 

ft 

ft 

ft 

0  3 

0  3 

•H 

O  i-a 

O  ft 

3 

S3 

S3 

S3 

3 

3 

3 

0 

0 

0 

0 

O 

ft 

ft 

N 

IS 

IS 

IS 

CO 

H  3 

3  ^ 

O 

O 

0 

O 

0 

3  0 

3  O 

3 

3 

3 

3 

ft 

CO  ft 

CO  ft 

ft 

ft 

ft 

ft 

0 

0 

0 

0 

£ 

3 

3 

ft 

ft 

ft 

£ 

£ 

£ 

£ 

ft 

3 

s 

0 

3 

3 

CO 

CM 

CO 

ft 

CO 

CO 

• 

• 

• 

• 

rH 

rH 

<M 

3 

CM 

On 

O 

<M 

NO 

S 

• 

• 

CO 

NO 

NO 

CM 

D- 

(H 

rH 

ft 

1 

1 

<0 

3 

O 

ft 

CO 

0 

ON 

On 

ON 

On 

1 — 1 

O- 

CM 

CM 

ON 

0 

• 

• 

• 

• 

• 

• 

C^ 

(M 

(M 

<M 

rH 

D 

1 — 1 

1 — 1 

CM 

(M 

CM 

CNJ 

4 

1 

1 

1 

1 

1 

1 

"0 

0 

1 — 1 

ft 

<3 

ft 

<3 

ft 

0 

ft 

£ 

1 

1 

1 

1 

1 

1 

3 

rn 

rH 

rH 

(M 

<M 

CM 

CM 

3  CD  CD 

3  o  o 

3 


£ 

0 

•H 

•H 

1 

O 

0 

0 

ft 

3 

> 

> 

ft 

3 

ft 

3 

3 

3 

0 

3 

0 

0 

0 

0 

0 

ft 

0 

ft 

ft 

ft 

ft 

ft 

ft 

3 

CO 

d 

CO 

3 

3 

ft 

S3 

rH 

S3 

rH 

3 

3 

3 

S3 

3 

•H 

3 

0 

3 

0 

ft 

ft 

0 

O 

rH 

O 

rH 

3 

3 

IS 

CO 

1 — 1 

•  H 

•H 

3 

3 

O 

S3 

•H 

ft 

0 

ft 

0 

O 

O 

3 

0 

O 

0 

O 

0 

O 

3 

H 

ft 

1 — 1 

CO 

3 

•H 

3 

•H 

hO 

bo 

£ 

£ 

£ 

£ 

s 

£ 

£ 

£ 

CM 

Cft 

-3- 

rH 

ON 

• 

• 

• 

• 

• 

• 

• 

• 

<T\ 

<r\ 

*0 

NO 

NO 

CO 

<M 

cj 

NO 

CM 

rH 

O 

0 

CM 

NO 

0- 

co 

CO 

1 — 1 

1 — 1 

rH 

1 — 1 

3 

0 

ft 

3 

3= 


0 

O 


CO 

CM 

CM 

-ft 

5 

8 

NO 

C^ 

O 

CM 

rH 

CO 

CO 

rH 

s 

• 

• 

• 

• 

• 

• 

• 

• 

CM 

i — 1 

rH 

1 — 1 

rH 

rH 

CM 

0 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

ft 

1 

1 

1 

1 

1 

1 

1 

I 

ft 

3 

3 

ft 

ft 

CO 


Cf 

< 

I 

p? 

ft 

1 

ft 

T 

< 

CM 

0^ 

<r^ 

r\ 

ft\ 

ft" 

ft" 

Table  5*  Oxygen  and  Hydrogen  Isotope  Ratios  in  Ground  Ice,  Surface  Water,  and  Groundwater 
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the  weighted  local  mean  precipitation,  and  have  not  undergone  isotopic 

TO 

enrichment  due  to  evaporation.  The  relatively  high  0  values  for  the 
surface  waters  suggest,  however,  that  high  summer  evaporation  rates  occur 
Examination  of  the  stable  isotope  profile  in  Figure  14  reveals 
three  unique  characteristics: 

’  18 

1.  the  frozen  peat  and  frozen  mineral  soil  0  values  differ  by 
approximately  1.0$#;  this  difference  is  significant  as  the 
analytical  technique  is  accurate  to  0.2$*. 

2.  within  each  group:  peat  and  mineral  soil;  stable  isotope  values 
are  nearly  constant;  values  for  reticulate  vein  ice  and  ice 
lenses  are  essentially  the  same. 

3.  the  isotopic  composition  of  the  sub-permafrost  groundwater  is 
approximately  1 .5%®  lighter  than  the  frozen  mineral  soil. 

It  is  unlikely  that  ground  ice  of  different  ages  is  accounting 
for  the  different  isotopic  ratios  in  the  peat  and  mineral  soil.  A 
previous  period  of  permafrost  degradation  is  not  indicated  by  the 
pal eo -temperature  determinations.  Oxygen  isotope  analysis  of  peat 
cellulose,  indicates  that  the  present  day  climate  is  the  warmest 
since  the  development  of  permafrost  in  this  region  (chapter  3). 

Further,  results  from  thaw  consolidation  tests  of  the  frozen  peat 
(chapter  5)  do  not  suggest  a  previous  thaw  period. 

The  difference  in  isotope  ratios  between  the  groundwater  and 
ice  may  be  due  to  fractionation  effects.  Fractionation  occurs 
during  the  freezing  of  water  under  equilibrium  conditions 
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(ie.  from  an  infinitely  large  source  of  well  mixed  water);  the  ice 

18 

containing  about  2  %o  more  0  than  the  water  from  which  it  has  frozen 
(Weston,  1955»  Amason,  1969)*  Mackay  and  Lavkulich  (1974)  and  Michels 
and  Fritz  (1978)  have  interpreted  isotope  differences  between  surface, 
subsurface  and  frozen  waters  for  an  aggrading  permafrost  site  as 
due  to  fractionation  effects. 

There  is  some  evidence  however  that  fractionation  may  not  have 
occurred  during  initial  freezing  of  the  paisa  soil,  but  rather  is  the  . 
result  of  post  freezing  modification  of  the  ice  core.  The  identical 
isotopic  values  for  reticulate  vein  ice  and  the  large  ice  lenses 
suggests  that  freezing  of  incoming  water  at  the  frost  interface 
occurred  very  rapidly  (Mackay, 1974;  McRoberts  and  Nixon,  1975)  with 
little  chance  for  isotopic  selection  to  occur.  Neither  is  a  Rayleigh 
fractionation  model  (Freidman  et  al,  1964)  supported  due  to  the  homogen¬ 
eous  isotopic  values  throughout  a  single  ice  lens.  A  second  possibility 
is  for  post  -  freezing,  "steady  state",  isotopic  exchange  to  have 
occurred  in  the  paisa.  Consider  an  ice  lens  system  with  unfrozen 
pore  water  beneath  the  lowest  ice  lens  (figure  16).  The  isotope  values 
of  the  frozen  soil  and  the  unfrozen  water  are  significantly  different, 
say  -30  %9  for  the  water  and  - 20 for  the  ice.  As  the  water  freezes 
onto  the  bottom  of  the  ice  lens,  water  will  also  melt  off-  the  top  of  the 
lens  due  to  temperature  induced  pressure  gradients  (Ershov,  1976; 

Mageau,  1978).  This  process,  called  regelation,  is  exactly  that 
outlined  by  Miller  (1972,1978)  for  pore  ice  movement  in  frozen  soils. 
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In  the  above  example,  the  water  freezing  at  the  bottom  of  the  lens 

will  become  isotopically  enriched  by  2  %o  to  become  -28  %*  .  Ice 

melting  at  the  top  of  the  lens  will  be  -20  %Q.  The  -28  %o  ice 

will  slowly  be  moved  up  the  ice  lens.  A  similar  process  involving  pore 

ice  regelation  and  an  unfrozen  water  flux  through  films  and  ice-free 

pores  (Miller,  1975)  will  occur  in  the  frozen  soil  above  the  lens. 

The  net  result  is  a  frozen  soil  with  near  homogeneous  isotope  profiles 

18 

in  the  ice,  underlain  by  unfrozen  water  which  is  depleted  2%0  in  0. 


In  summary ,  stable  isotopic  analysis  of  the  paisa  cores  has  shown 
that  there  was  minimal  evaporation  of  the  soil  water  prior  to  its 
being  permafrozen.  This  supports  the  conclusion,  made  in  a  later  chapter, 
that  pal eo temperatures  were  significantly  cooler  than  present  day  temper¬ 
atures.  Separate  groundwater  flow  systems  for  the  peat  and  mineral  soils 
were  earlier  deduced  from  the  water  chemistry.  This  conclusion  was 
supported  by  the  isotope  analysis.  Tritium  dating  revealed  an  anon- 
omously  young  age  for  the  ground  ice  as  compared  to  the  paisa  itself,  . 
suggesting  some  water  exchange  has  taken  place.  This  conclusion  was 
supported  by  analysis  of  the  stable  isotope  profiles,  where  a  model 
was  proposed  for  regelation  flow,  with  accompanying  isotope  exchange, 
as  water  moves  upward  in  response  to  thermal  gradients. 


' 
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2.10  Hydrogeology  of  the  Paisa  Complex 

An  investigation  was  made  of  the  groundwater  regime  within  the 
paisa  complex  for  two  reasons.  First,  to  document  the  pore  pressures 
immediately  beneath  a  paisa.  Secondly,  to  determine  whether  the  col¬ 
lapse  bogs  and  thermokarst  sinkholes  are  underlain  by  permafrost  at 
depth  or  are  local  sources  of  groundwater  recharge.  Recharge  may  be 
expected  as  water  levels  in  the  bogs  are  higher  than  the  surrounding 
fen  (see  Figure  4)  with  no  apparent  surface  drainage. 

For  each  collapse  feature  studied,  piezometers  were  installed  by 
driving  well  points  in  the  deepest  part  of  the  bog  or  pond  in  order  to 
eliminate  two-dimensional  effects.  Piezometers  that  pierced  the  perma¬ 
frost  were  installed  during  the  winter  drilling  program.  The  saturated 
bogs  acted  as  a  body  of  standing  water  with  no  head  loss  with  depth. 
Piezometer  installations  were  completed  in  the  usual  manner  then 
charged  with  a  mixture  of  kerosene  and  methyl  alcohol  (unit  weight 
equal  to  water)  in  order  to  prevent  freezing.  Hydraulic  conductivities 
were  determined  by  the  Hvorslev  (1949)  water  level  response  technique, 
using  the  'University  of  Alberta  transducer  piezometer'  (see  McRoberts 
and  Morgenstem,  1974  for  full  description).  Head  measurements  are  re¬ 
ported  relative  to  the  level  of  the  water  surface  of  the  surrounding 
fen. 


As  would  be  expected  in  a  discontinuously  frozen  area,  no  excess 
hydrostatic  heads  were  recorded.  All  piezometric  measurements  were 
consistent  in  showing  a  decrease  in  fluid  potential  with  depth.  The 
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head  measurements  and  approximate  isopotentials  for  one  of  the  collapse 
hogs  is  illustrated  in  Figure  17  .  The  decrease  with  depth  indicates 
that  the  collapse  features  are  not  underlain  hy  permafrost ;  that  water 
is  moving  from  the  raised  hog,  beneath  the  paisa,  and  discharging  into 
the  surrounding  fen. 

Estimates  of  the  volume  of  groundwater  supplied  hy  the  raised 
hogs  can  he  determined  from  Darcy's  equation: 

Q  =  k  i  A 

where  Q  is  the  volume  of  flow 

—7 

k  is  the  permeability,  which  ranged  from  2.9  x  10  '  cm/sec  to 

—7  —7 

6.8  x  10  cm/sec;  averaging  5»5  x  10  cm/sec 

i  is  the  hydraulic  gradient;  averaging  .08 

2 

A  is  the  hog  area;  300  m 

The  calculated  recharge  of  4  ittV  year  corresponds  to  a  water 
thickness  of  1.3  cm  over  the  hog  surface  or  roughly  5  %  of  the  mean 
annual  precipitation.  Assuming  a  steady  water  table,  evaporation  must 
nearly  equal  precipitation  for  the  closed  hogs.  Oxygen  isotope  anal¬ 
ysis  of  hog  surface  waters  (Chapter  4  )  confirms  the  high  evaporation 
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FIGURE  1 V  ;  ISOPOTENTIALS  AND  LOCAL  GROUNDWATER  FLOW  DUE 

PERCHED  SATURATED  COLLAPSE  BOG;  SITE  E-3  . 


CHAPTER  3:  PALEOCLIMATIC  DETERMINATION 
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3.1  Introduction 


Paleoclimatic  information  is  vital  in  the  determination  of  the 
history  of  permafrost  aggradation  and  degradation.  While  other 
elements  of  a  regional  climate  play  a  role  in  permafrost  development, 
temperature  is  acknowledged  as  the  most  important  component  (Brown, 
1966) .  In  the  Arctic  and  Sub-Arctic  regions  of  North  America,  paleo- 
temperatures  have  been  estimated  from  palynological  studies,  peat 
stratigraphy,  and  freshwater  paleontology.  All  of  these  studies  have 
relied  on  biological  response  to  inferred  climatic  change.  The  de¬ 
pendence  of  biological  systems  on  many  local  edaphic  factors  apart 

from  climate  is  well  documented  (Iverson,  1953)  the  result  being  that 
many  of  these  studies  have  been  contradictory  and  yielded  little  quant¬ 
itative  data. 

The  determination  of  oxygen  and  deuterium  isotope  ratios  in  non¬ 
exchangeable  plant  cellulose  has  been  recently  shown  to  reflect  the 
isotope  ratio  of  the  water  used  by  the  plant  in  its  synthesis  of  cel¬ 
lulose  (Epstein  et  al,  1976;  Thompson  and  Gray,  1977).  Numerous 

studies  (Epstein  and  Mayeda,  1953;  Dansgaard,  1964;  Friedman,  1964) 

l8  l6 

have  shown  that  d/h  and  iO0/  °0  ratios  in  meteoric  waters  reflect 
the  temperatures  at  the  sites  of  precipitation;  the  lower  the  temp¬ 
erature,  the  smaller  the  isotopic  ratio.  The  climatic  significance  of 

18  l6 

d/r  ratios  (Epstein  and  Yapp,  1976)  and  1  0/  “  0  ratios  (Gray  and 
Thompson,  1977)  have  subsequently  been  demonstrated  for  tree-ring 
sequences  from  different  climatic  regimes.  It  follows  that  the  iso¬ 
topic  variations  in  plant  cellulose  may  reflect  the  temperature  of 


. 
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precipitation  at  the  growth  site. 


In  this  chapter,  oxygen  isotope  ratios  of  peat  cellulose  are  used 
to  estimate  paleotemperatures  at  the  study  site  over  the  last  10,000 
years.  This  information  will  be  used  in  later  discussion  of  the  his¬ 
tory  of  permafrost  in  Upper  Mackenzie  Valley  peatlands. 


3.2  Previous  Investigations  of  Holocene  Paleoclimate  in  N.  W,  North 

America 

A  comprehensive  review  of  paleoclimatic  investigations  in  the 
Arctic  and  Sub-Arctic  regions  of  North  America  has  been  compiled  by 
Ritchie  and  Hare  (l97l).  The  reviewers  found  many  of  the  results, 
especially  that  based  on  pollen  evidence,  to  be  contradictory.  While 
some  workers  supported  a  mid-Holocene  warm  period  (Livingstone,  1957 ; 
Heusser,  i960);  others  have  concluded  that  "the  climate  has  never 
been  significantly  warmer  than  it  is  at  the  present"  (Colinvaux,  1967) . 
Megafossil  evidence  was  more  definitive.  Definite  evidence  of  the 
tundra-taiga  boundary  being  significantly  further  north  than  at  present, 
suggesting  greater  mean  summer  temperatures,  has  been  found  in  the  form 
of  buried  wood  and  fossil  beaver  artifacts  in  the  following  areas: 

Seward  Peninsula  (McCulloch  and  Hopkins,  1966);  Kendal  Island  (Mackay, 
1971);  Tuktoyaktuk  Peninsula  (Ritchie  and  Hare,  1971);  and  in  buried 
forest  paleosols  at  Keewatin  (Bryson  et  al,  1965) .  Notable  studies 
carried  out  in  recent  years  include  new  results  by  Ritchie  and  Hare 
(1971);  Matthews  (1974);  Rampton  (l97l);  Zoltai  and  Tarnocai  (1975); 
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and  Delorme  et  al  ( 1977 ) . 


The  pollen  stratigraphy  of  Kettle  Lakes  near  Tuktoyaktuk  is 
interpreted  by  Ritchie  and  Hare  (l97l)  as  indicative  of  mean  summer 
termperatures  approximately  5  °  C.  higher  than  at  present,  with  the 
growing  season  about  30  days  longer,  during  the  time  interval  8 500 
years  to  5500  years  Before  Present.  J.R.  Mackay  (personal  commun¬ 
ication)  feels  that  the  large  summer  temperature  increase  at  their 
site,  may  be  more  apparent  than  real.  During  this  time  interval  sea 
level  was  at  least  15  m  below  present  day  levels  (Mackay,  1963) .  Bath¬ 
ymetric  contours  off  the  Tuktoyaktuk  peninsula  indicate  that  a  rel¬ 
ative  sea  level  lowering  of  15  m  would  result  in  a  coast  line  dis¬ 
placed  50  km  to  the  north.  By  assuming  temperature  gradients  similar 
to  those  existing  today,  much  of  the  mean  summer  temperature  increase 
postulated  by  Ritchie  and  Hare  could  be  accounted  for  by  the  increased 
distance  from  the  Arctic  Ocean. 

Matthews  (l97^)»  in  a  palynological  study  of  colluviated  silts 
near  Fairbanks,  Alaska  found  evidence  for  arid  steppe-tundra  conditions 
from  approximately  12,000  years  B.P.  until  sometime  between  9200  and 
7800  years  B.P.  He  postulated  that  Central  Alaska  possessed  an  arctic 
climate  similar  to  present  day  Northern  Alaska  during  this  time  in¬ 
terval.  Alder  fluctuations  in  the  pollen  record  suggest  that  climatic 
warming  followed,  leading  to  a  taiga  vegetation. 

Rampton  (1971),  has  attempted  to  quantify  his  pollen  record  from 
the  Snag-Klutlan  area,  Yukon  Territory.  The  former  vegetative 
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community  suggests  “both  mean  July  temperatures  and  mean  annual  temp¬ 
eratures  approximately  7°  to  8°  cooler  than  present  values  during  the 
time  interval  13,500  to  10,000  years  B.P.  Mean  July  temperatures  and 
mean  annual  temperatures  of  approximately  4.5°  C.  to  5*5°  0.  cooler 
than  present  values  were  proposed  for  the  time  interval  10,000  years 
to  8700  years  B.P.  The  pollen  evidence  indicates  similar  vegetation 
to  the  present  day  since  8700  years  B.P.  The  presence  of  fossil  logs 

above  tree  line  suggest  that  temperatures  may  have  fluctuated  up  to 
1°  C.  above  present  day  values  between  6000  years  and  1220  years  B.P. 

Tree  ring  studies  indicate  temperatures  during  the  last  200  years  were 

at  least  1°  cooler  than  present. 

Zoltai  and  Tarnocai  (1975)  and  Zoltai  and  Pettapiece  (1978)  have 
used  peat -permafrost  relationships  and  dated  organic  horizons  in 
cryoturbated  soils  to  suggest  that  during  the  time  interval  8000  years 
to  4000  years  B.P.,  permafrost  was  much  less  widespread  throughout  the 
Mackenzie  Valley.  They  believe  that  a  climatic  cooling  some  3000  to 
4000  years  B.P.  initiated  permafrost  conditions,  restricting  peat 
accumulation  to  unfrozen  fens. 

Most  recently,  Delorme  et  al  (1977)  have  used  ostracod  assemblages 
to  infer  climatic  conditions  in  the  Mackenzie  Delta  -  Northwestern 
Arctic  Coast  area  during  the  time  interval  14,410  to  6820  years  B.P. 
Mean  annual  temperatures  of  up  to  +  1.2°  C.,  approximately  10  G. 
higher  than  present  day  values,  and  annual  precipitation  from  55  to 
234  mm  greater  than  the  present  day,  were  postulated. 
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Mackay  (1978),  in  response  to  the  paper  by  Delorme  et  al,  has 
cited  both  the  preservation  of  icy  Pleistocene  sediments  and  the  lower 
portions  of  relic  ice  wedges  beneath  a  widespread  thaw  unconformity 
along  the  Arctic  Coast  as  evidence  that  maximum  postglacial  thaw  was 
only  a  few  meters.  The  depth  of  thaw  indicates  some  post-glacial 
climatic  warming,  however,  mean  annual  temperature  must  have  remained 
several  degrees  below  0°  C.  during  this  time  interval. 

In  summary,  there  is  reasonably  strong  evidence  for  a  post-glacial 
period  in  North  Western  Arctic  Canada  when  tree  lines  were  further 
north  and  the  active  layer  was  thicker  than  at  the  present  time.  These 
changes  were  no  doubt  due  to  an  increase  in  mean  summer  temperatures. 
Temperature  increases  of  the  magnitude  suggested  by  Ritchie  and  Hare, 
as  well  as  Delorme  et  al  are  not  substantiated.  Evidence  for  mean 
summer  temperatures  substantially  greater  than  those  of  the  present 
day  is  not  as  compelling  for  inland  sties  as  it  is  for  arctic  coastal 
sites. 

3.3  Methods  of  Analysis 

Frozen  peat  samples  were  collected  with  a  CRREL  (Cold  Regions 
Research  and  Engineering  Laboratory  -  U.S.  Army)  barrel  coring  device 
form  site  E  11  within  the  frozen  paisa  complex. 

Cellulose  was  extracted  from  the  peat  samples  using  a  modification 
of  a  method  described  by  Theander  ( 195*0 »  which  is  illustrated  by  a 
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flow  chart  in  Figure  18  .  The  method  yields  a  greyish-white  to  light 
grey  substance.  An  X-ray  diffraction  pattern  of  the  substance  reveals 
strong  peaks  in  the  range  characteristic  of  alpha-cellulose.  Following 
extraction,  approximately  25  mg  of  vacuum  dried  sample  is  weighed  into 
a  nickel  boat,  placed  under  vacuum,  then  combusted  in  a  reaction  ves¬ 
sel  at  1150  0  G.  for  1  hour  (see  Thompson  and  Gray,  1976).  The  reac¬ 
tion  produces  a  mixture  of  CO,  CO^  and  gasses.  The  gas  dif¬ 
fuses  through  the  walls  of  the  reaction  vessels,  forcing  the  reaction 
to  completion.  After  combustion,  the  gasses  are  trapped  in  a  silica 
gel  trap  maintained  at  liquid  nitrogen  temperature.  The  GO  fraction 
is  converted  to  GO^  and  G  by  sparking  between  platinum  electrodes.  The 
GO^  fractions  are  combined  and  the  yields  measured  mano metrically.  The 
peat  samples  averaged  90  to  of  the  theoretical  predictions,  indi¬ 
cating  trace  amount  of  contaminant.  All  peat  samples  were  run  as  rep¬ 
licates;  if  replicate  results  did  not  agree,  a  third  sample  was  run. 
The  precision  of  the  replicate  analysis  was  better  than  -  Iso¬ 

tope  measurements  were  made  on  a  90  °  sector,  25  cm  double  collecting 
isotope  ratio  mass  spectrometer. 

3.4  Results 

Isotope  measurements  are  expressed  as  a  ^  value  (°/oo)  with  re¬ 
spect  to  Standard  Mean  Ocean  Water  (SMOW).  The  results  are  listed  in 
Table  6  . 


The  isotopic  results  are  plotted  with  respect  to  depth  in  Figure 
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FIGURE  18 


FLOW  CHART  FOR  PEAT  CELLULOSE  EXTRACTION 


Peat  (air  dry) 


( soluble) 


X 


Waxes,  resin, 
tannins,  bitumens, 
hydrocarbons 


( soluble) 


Benzene/Methanol  (2:l)  Extraction, 
6  hours 


Hot  water  (?03  C)  3*3  hours 


-,0  p 
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TABLE  6 

PEAT  OXYGEN  ISOTOPE  RESULTS 
-  Core  E  3 


Sample 

Depth  (m) 

Analysis  number 

180 

SMOW 

Standard 

deviation 

E  3-1 

.05 

3 

'  16.80 

.40 

E  3-2 

.20 

2 

17.02 

.14 

E  3-3 

.51 

2 

16.58 

.39 

E  3-4 

.77 

2 

16.17 

.10 

E  3-3 

1.00 

2 

16.66 

.33 

E  3-6 

1.25 

2 

16.44 

.32 

E  3-7 

1.42 

2 

16.05 

.19 

E  3-8 

1.63 

2 

15.38 

.40 

E  3-9 

1.82 

3 

15.85 

.15 

E  3-10 

2.12 

2 

16.36 

.35 

E  3-11 

2.30 

3 

16.84 

.20 

E  3-12 

2.50 

2 

16.45 

.31 

E  3-13 

2.75 

4 

15.55 

.2? 

E  3-14 

3.15 

2 

15.70 

.15 

E  3-15 

3M 

3 

15.50 

.41 

' 

19  ,  together  with  a  radio -carbon  age  depth  profile  and  an  illustra¬ 
tion  of  the  peat  stratigraphy.  A  description  of  the  peat  character¬ 
istics  and  phytological  components  can  be  found  in  chapter  1.  The 
radiocarbon  plot  depicts  corrected  C-l4  ages  as  determined  by  the 
University  of  Waterloo  Isotope  Lab  (Wat  400  to  405).  In  constructing 
the  age-depth  profile,  rates  of  peat  accumulation  were  assumed  linear 
between  dated  samples.  Peat  accumulation  began  10,380  years  B.P., 
terminating  approximately  400  years  B.P.  Analysis  of  modem  peat 
forming  plants  at  the  study  site  has  allowed  extrapolation  of  the  iso¬ 
tope  curve  to  the  present  day. 

A  possible  hiatus  due  to  a  period  of  freezing  then  uplift  or  to 

drying  of  the  peatland  is  found  at  a  depth  of  1.6  m.  This  level  is 

\ 

marked  by  a  2  cm  layer  of  woody  nanolignid  peat,  similar  to  the  forest 
peat  layer  found  on  the  surface  of  modern  paisas.  The  interval  was 
apparently  short  lived;  a  plot  of  the  C-l4  age  of  the  woody  peat  lies 
only  slightly  off  the  line  connecting  peat  ages  on  either  side  of  the 
layer  ( Figure  19  ) . 

3.41  Effect  of  Different  Plant  Species 

Prior  to  analysing  the  isotopic  composition  of  bulk  peat 
samples,  it  was  necessary  to  determine  whether  individual  peat-forming 
plant  species  have  similar  oxygen  isotope  ratios.  Samples  of  living 
plants,  representing  all  the  major  peat  forming  species,  were  col¬ 
lected  from  the  study  site  during  the  summer  field  season.  The 
oxygen  isotopic  composition  of  each  of  the  species  was  determined;  the 
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results  are  listed  in  Table  7  . 


TABLE  7 

Isotopic  Composition  of  Modern  Peat  Forming  Plants 


Plant  Species 

180 

Sphagnum  fuscum 

18.65 

Sphagnum  recurvum 

18.51 

Drepanocladus 

18.37 

Carex  aquatilis 

18.40 

Carex  (various  spp.) 

18.47 

Duckweed 

18.39 

Mixture 

18.48 

SMOW 

“  0.2 

%o 

%o 

%o 

%0 

%o 

%o 

%o 


The  results  show  that,  within  the  limits  of  analytical  error,  all 
the  analyzed  plant  species  have  very  similar  oxygen  isotopic  values, 
reflecting  the  constant  isotopic  composition  of  the  water  in  which 
they  grew.  It  is  therefore  not  necessary  to  separate  individual 
species  from  the  peat;  bulk  peat  samples  will  give  representative 


results. 


3.42  Effect  of  Humification  on  Isotope  Values 


There  is  a  possibility  of  oxygen  exchange  between  the  C-0 
bonds  in  the  peat  cellulose  and  percolating  bog  waters  or  with  bacteria 
during  decomposition.  If  such  exchange  did  occur,  the  result  would  be 
to  produce  an  homogenization  of  isotope  values  in  the  peat.  Evidence 
suggests,  however,  that  such  isotope  exchange  is  very  unlikely;  the 
C-0  molecular  bond  in  cellulose  being  very  strong. 

No  oxygen  exchange  was  detected  in  cellulose  samples  that  were 
subjected  to  harsh  oxidizing  conditions  (Gray  and  Thompson,  1976; 
Epstein,  Yapp,  and  Hall,  1976;  this  thesis)  or  to  high  partial  pres¬ 
sures  of  oxygen  (this  thesis).  Further,  Yapp  and  Epstein  (1977) 
found  no  evidence  for  isotope  exchange  having  occurred  in  wood  that 
had  been  buried  in  wet  soils  for  11,800  years.  In  conclusion,  the 
evidence  to  date,  does  not  support  the  probability  of  isotope  exchange 
occurring  in  the  peat  during  humification. 

3.43  Temperature  Effects 

The  effect  of  temperature  on  the  ultimate  isotope  values  of 

plant  cellulose  is  very  complex.  The  temperature  of  air  mass  vapor 

condensation  largely  determines  the  isotopic  composition  at  a  specific 

site,  but  the  degree  of  proportionality  is  a  function  of  the  history 

of  the  air  mass  (Dansgaard,  1964;  Friedman  et  al,  1964).  The  effect 

of  evaporation  and  evapo transpiration,  both  temperature  dependent,  on 
1  ft 

increasing  0  values  in  plant  cellulose  has  been  discussed  by  Epstein 


« 


et  al  (1977).  When  comparing  localities  that  have  widely  different 

humidities,  this  effect  will  he  of  importance.  Gray  and  Thompson 

(1976)  however  have  shown  that  at  a  specific  locality,  where  large 

humidity  changes  over  time  are  not  usually  encountered,  excellent 

l8 

agreement  can  be  found  between  Mean  Annual  Temperature  and  0  (cellu¬ 
lose).  In  peat  bogs,  humidity  is  always  high,  and  effects  due  to 
changes  of  evapo transpiration  rates  can  probably  be  neglected,  es¬ 
pecially  when  examining  variations  at  a  specific  site.  If  enrich¬ 
ment  due  to  evapotranspiration  has  occurred,  the  effect  will  be  to 
overestimate  Mean  Annual  Temperatures.  Fractionation  occurring 
during  cellulose  synthesis  may  also  be  temperature  dependent  (Epstein 
et  al,  1977;  Gray  and  Thompson,  1976).  In  all  previous  site  specific 
studies,  the  temperature  coefficient  for  cellulose  was  positive,  ranging 
between  1  and  3  °/oo  ^0/  °G.  (Libby  and  Pandalfi,  197^;  Epstein  and 
Yapp,  1977;  Gray  and  Thompson,  1976). 

3.5  Paleoclimatic  Interpretations 

The  preceding  discussion  documented  that  Mean  Annual  Air  Temp¬ 
eratures  are  strongly  correlated  with  oxygen  isotope  ratios  in  plant 
cellulose.  Before  pal eo -temperature  interpretations  can  be  made  of  the 
isotope  fluctuations  however,  a  temperature  coefficient  must  be  found 
for  peat  cellulose.  Unfortunately,  peat  cannot  be  calibrated  in  the 
same  way  as  tree  rings  because  it  is  impossible  to  separate  yearly 
accumulations.  The  temperature  coefficient  for  Edmonton  Spruce,  1.3°/oo 
^0/  0  G.  (Thompson  and  Gray,  1976) ,  mi^it  be  applicable;  however 
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Sphagnum  cellulose  may  fractionate  differently  than  spruce. 


A  method  of  comparing  Sphagnum  and  tree  ring  temperature  coef- 

18 

ficients  is  to  plot  ^(cellulose)  ^an^s  ^rom  various  locations 
versus  the  Mean  Annual  Temperature  at  that  site.  Sphagnum  moss  from 
a  wide  variety  of  climatic  zones  was  analyzed,  then  results  plotted 
against  Mean  Annual  Temperature  (average  of  last  5  years)  in  Figure  20. 
For  some  of  the  locations,  Sphagnum  samples  were  available  from  more 
than  one  bog;  the  range  of  results  at  these  locations  is  plotted  as 
a  bar.  On  the  same  graph,  oxygen  isotope  values  of  modern  tree  ring 
cellulose  from  various  sites  are  plotted  (unpublished  data,  P. 

Thompson) .  This  data  set  includes  the  previously  mentioned  Edmonton 
spruce  analyzed  by  Gray  and  Thompson  (1976).  The  isotopic  values  from 
northern  continental  North  America  (Edmonton,  Alberta;  North  Bay, 
Ontario;  Fort  Simpson,  N.W.T.;  Old  Grow,  Yukon;  and  Coppermine, 
N.W.T.)  are  reasonably  defined  by  the  linear  relationship: 

<$l80  =0.5T  +  19.90 

The  values  for  subtropical  and  tropical  localities  (Houston,  Texas; 
Gampinos,  Brazil)  fall  significantly  below  this  line,  while  values 
from  a  Pacific  maritime  area  (Queen  Charlotte  Islands,  B.C.)  fall  above 
the  line. 

The  tree  ring  isotope  values  respond  in  a  similar  manner;  the 
regression  line  is  parallel  to  the  peat  line  but  displaced  2  °/oo 
higher  (Figure  20  ).  The  reason  for  the  displacement  is  not  apparent 
but  may  be  due  to  a  degree  of  isotopic  equilibration  of  atmospheric 
C09  with  the  tree  water  prior  to  photosynthesis  and  the  fixation  of 
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Figure  20  :  Change  in  oxygen  isotope  content  of  modern  Sphagnum  moss  with  the  mean 
annual  temperature  of  their  place  of  growth. 


COp  into  cellulose  (Gray  and  Thompson,  1976). 


In  theory,  this  method  should  lead  to  a  temperature  coefficient 
for  plant  cellulose.  Comparison  of  Gray  and  Thompsons'  tree  ring 
coefficient  of  1.3  °/oo/  °C.  with  the  value  of  0.5  °/oo  /  0  C.  from 
Figure  20  ,  clearly  shows  this  is  not  the  case.  As  previously  dis¬ 
cussed,  many  non-temperature  related  factors  can  effect  isotopic 
fractionation.  These  effects  become  apparent  when  comparing  different 
localities,  hence  also  the  scatter  of  the  data  points,  especially  for 
more  tropical  samples.  The  graph  does  show  that  peat  and  tree  ring 
isotope  ratios  increase  in  parallel  as  Mean  Annual  Temperatures  in¬ 
crease.  While  the  need  for  a  great  deal  more  work  is  acknowledged,  use 
of  the  temperature  coefficient  for  White  Spruce  (1.3  °/°°  /  °C.)  should 
give  the  approximate  magnitude  of  temperature  fluctuations  during  peat 
accumulation. 

The  oxygen  isotope  curve,  from  analysis  of  peat  core  E  11  (Figure 
19),  is  characterized  by  relatively  little  variation  throughout  the 
profile.  A  comparatively  cool  climate  is  indicated  between  approx¬ 
imately  10,000  years  B.P.  and  7000  years  B.P.  This  was  followed  by  a 
warming  trend  between  approximately  6000  years  B.P.  and  3000  years  B.P. 
While  this  time  interval  corresponds  to  that  of  the  classical  Hypsi- 
thermal,  the  isotope  values  do  not  indicate  exceptional  warming.  A 
cooling  trend  followed,  reaching  a  minimum  value  between  approximately 
4000  and  3000  years  B.P.;  this  was  the  lowest  isotope  value  measured 
in  the  entire  peat  column.  Interestingly,  a  change  in  peat  type  is 
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encountered  immediately  above  this  layer,  from  a  Sphagnum  peat  to  a 
forest  peat,  suggesting  a  brief  period  of  permafrost  invasion  then 
colonization  by  Black  Spruce.  Isotope  variations  are  small  in  the 
upper  1.5  meters  of  core;  a  gradual  warming  trend  is  indicated, 
especially  over  the  past  400  years. 

Modern  peat  producing  plants  in  the  Fort  Simpson  area  have  higher 
isotope  ratios  than  any  of  the  peat  samples.  This  suggests  current 
Mean  Annual  Temperatures  are  the  warmest  that  have  occurred  in  the  past 
10,000  years.  Using  the  above  temperature  coefficient  of  1.3  °/oo/  °C., 
a  minimum  Mean  Annual  Temperature  approximately  3°  0.  colder  than 
present,  or  -7°  0.,  is  indicated. 

An  independent  method  of  checking  these  temperature  approximations 
is  to  examine  nearby  deep  temperature  profiles.  At  depth,  the  geo¬ 
thermal  gradient  is  slow  to  respond  to  surface  temperature  changes; 
hence  perturbations  in  the  gradient  can  be  used  to  calculate  previous 
ground  surface  temperatures  (Judge,  197*0  •  A  deep  temperature  profile 
from  an  oil  exploratory  well,  200  km  east  of  Fort  Simpson,  is  plotted 
in  Figure  21  (data  from  Judge,  197*0.  The  well  penetrates  33*5  m  of 
unfrozen  unconsolidated  material,  underlain  by  300  m  of  shale  and  113  m 
of  carbonates.  The  temperature  gradient  is  abnormally  low  to  a  depth 
of  66  m,  which  does  not  correspond  to  any  change  in  lithology,  hence 
is  indicative  of  recent  climatic  warming. 


Extrapolation  of  the  lower  unaffected  part  of  the  curve  to  the 
surface  indicates  a  Paleo-Mean  Annual  Surface  Temperature  (M.A.S.T.)  of 


■  , 


DEPTH  (meters) 


76 


77 


-2.8  °G.  Compared  to  the  present  day  M.A.S.T.  of  +  0.4  °C.,  this 
represents  a  warming  of  +3*2  °C.  Prior  to  the  onset  of  surface 
warming,  at  least  30  m.  of  permafrost  was  present  (Figure  21).  Degrad¬ 
ation  of  the  permafrost  was  by  bottom  melting  due  to  the  geothermal 
heat  flux.  The  rate  of  melting  of  the  permafrost  due  to  the  geo¬ 
thermal  heat  flux  is  calculated  to  be  1.0  to  1.5  cm/  year,  hence 
surface  warming  must  have  commenced  roughly  2000  to  3000  years  B.P. 

A  Mean  Annual  Temperature  approximately  3  °C.  cooler,  2500  years  B.P., 
agrees  well  with  the  paleotemperature  predictions  based  on  oxygen  iso¬ 
topes  from  the  peat  cellulose. 
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CHAPTER  4:  EVOLUTION  OF  THE  PEATLAND  LANDSCAPE 
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4.1  Introduction 


The  evolution  of  a  peatland  is  recorded  in  the  peat  stratigraphy. 
The  interrelationships  of  mineral  soil  topography,  hydrologic  regimes, 
vegetational  sequences,  hog  processes  and  climate  create  distinct 
peatland  types;  therefore,  the  peat  is  more  than  just  a  record  of 
past  plant  communities.  The  work  of  Kulczynski  (1949)  and  Heinselman 
(1963)  is  unique  in  examining  past  environmental  conditions,  as  de¬ 
duced  from  the  peat  stratigraphy,  for  their  role  in  creating  the  mod¬ 
ern  peatland  landscape.  This  type  of  integrated  study  has  not  "been 
carried  out  in  discontinuously  frozen  peatlands. 

In  this  chapter,  the  change  in  physical  environmental  conditions 
in  the  Fort  Simpson  peatland  over  the  past  10,000  years  is  inferred 
from  the  peat  stratigraphy.  The  stratigraphy  across  the  paisa  complex 
records  six  vegetative  communities: 

1.  organic  detritus  in  silt 

2 .  aquatic  plants 

3.  brown  moss-sedge 

4.  sedge-shrub-Black  Spruce 

5.  Sphagnum  bog 

6.  Black  Spruce-Cladonia-feathermoss 

Each  of  these  communities  implies  a  different  wetland  environment. 
Modern  analogues  of  these  communities  allow  reconstruction  of  en¬ 
vironmental  conditions  at  the  time  of  deposition.  Six  radiocarbon 
dates  permit  a  chronological  sequence  to  be  outlined  for  the  vegetative 
and  environmental  successions. 


. 


4.2  Environmental  Sequences 


Immediately  after  the  retreat  of  the  Laurentide  Ice  Sheet,  Glacial 
Lake  McConnell  occupied  large  areas  of  the  Mackenzie  Valley  (Rutter, 
1974) .  While  the  lake  was  subject  to  numerous  fluctuations  creating 
complex  beach  ridges  around  its  margins,  deep  water  conditions  per¬ 
sisted  until  approximately  11,000  years  B.P.  Pollen  profiles  indicate 
that  during  the  approximately  1000  year  existence  of  the  lake,  boreal 
forest  vegetation  occupied  the  poorly  drained,  possibly  discontin- 
uously  frozen,  glacial  soils  bordering  the  lake  (T.  Habgood,  personal 
communication).  Black  Spruce  was  the  dominant  tree  species,  being 
probably  the  first  to  re-colonize  the  area  as  it  is  best  adapted  to 
poorly  drained  soils.  Organic  deposition  in  the  turbid  waters  was  very 
slow  but  created  a  relatively  thick  layer  of  organic  detritus  in  a  silt 
matrix. 

Drainage  and  evaporation  of  Lake  McConnell  was  not  sudden  but  in¬ 
volved  a  progressive  lowering  of  base  levels.  Eventually,  a  level 
plain  with  many  shallow  lakes  and  ponds  remained.  This  phase  was  rel¬ 
atively  short  lived  and  resulted  in  the  discontinuous  deposition  of 
shallow  organic  gyttja  and  aquatic  peat. 

It  is  uncertain  whether  the  lake  plain  completely  dried  up.  A 
dessication  period  has  been  reported  in  northern  Alberta  (May  and 
Thomson,  1977)  for  this  time  interval.  Further,  near  surface  lacus¬ 
trine  sediments  are  overconsolidated  (see  section  3 A),  a  character¬ 
istic  of  dessicated  soils  (Wu,  1977).  Also,  calcium  sulphate  concen- 
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trations  were  recorded  at  the  top  of  the  mineral  soil  (Chapter  l),  a 
characteristic  of  evaporated  lake  waters.  Conflicting  evidence  is  the 
lack  of  a  phyto  response  in  the  pollen  diagrams  from  the  area.  If  an 
arid  interval  coincided  with  the  drainage  of  Lake  McConnell,  it  was 
probably  of  short  duration. 

Fens  quickly  overspread  the  nutrient  rich  lacustrine  plain  and 
paludification  of  the  landscape  began.  Fen  communities  can  form  peat 
capable  of  holding  back  water  hence  causing  significant  rises  in  the 
water  table  (Heinselman,  1970  ) .  This  event  is  dated  by  the  basal  layer 
of  brown  moss  and  sedge  fen  peats  that  nearly  uniformly  cover  the  form¬ 
er  glacial  lake  basin,  dated  at  approximately  8 500  years  B.P.  Radio¬ 
carbon  dated  horizons  and  the  uniform  thickness  of  the  brown  moss- 
sedge  peat  layer  suggests  that  this  community  lasted  for  several  thou¬ 
sand  years.  The  area  around  the  large  lakes  to  the  north  of  the  paisa 
complex  is  a  modern  analogue  of  the  treeless,  very  gently  sloping, 
minera trophic  fens  that  occupied  the  site. 

As  the  peat  accumulated,  drier  sites  became  available.  Shrubs 
such  as  Sweet  Gale,  dwarf  birch,  willow  and  to  a  lesser  extent  trees 
like  Black  Spruce  and  Tamarack  invaded  the  fens;  in  some  localities 
a  shrub  fen  or  carr  resulted.  The  vegetative  community  was  floristic- 
ally  rich:  the  peat  contains  sedges,  brown  mosses,  shrubs,  wood, 
mosses  and  possibly  grasses.  Spruce  is  much  more  abundant  in  the 
woody  fen  peat  than  that  described  by  Zoltai  and  Tarnocai  (1975)  from 
peat  plateaus  in  the  central  Mackenzie  Valley. 


The  thick  Sphagnum  peat  layer  found  in  all  cores  is  testimony  to 
a  change  from  nutrient  rich  mineratrophic  fens  to  nutrient  poor  ombro- 
trophic  bogs.  Prior  to  the  Sphagnum  bloom,  the  fen  surface  had  a  very 
gentle  uniform  slope;  the  continual  accumulation  of  peat  tending  to 
level  out  the  peatland  surface.  Ombrotrophic  conditions  developed 
when  the  peatland  surface  rose  above  the  level  of  percolating,  mineral 
rich,  groundwater;  hence  creating  a  local  water  table  divide.  The 
divide  may  have  been  partly  present  in  the  original  physiography,  or 
was  developed  entirely  by  peat  accumulation,  or  was  created  by  stream 
capture  (Heinselman,  1970) .  The  importance  of  local  water  table  div¬ 
ides  in  the  development  of  large  Sphagnum  bogs  is  also  stressed  by 
Kulc  zynski  ( 1949 ) . 

The  lowering  of  regional  water  tables,  possibly  related  to  the 
more  integrated  drainage  system  that  was  developing,  initiated  peat¬ 
land  drainage  in  more  than  one  direction.  Areas  of  elevated  mineral 
soil  beneath  the  peat  provided  the  first  areas  to  develop  local  water 
table  divides.  The  lake  marginal  area  was  particularly  favored,  hav¬ 
ing  an  irregular  substratum  topography  (Chapter  l).  Sphagnum  quickly 
invaded  these  ombrotrophic  sites  creating  extensive  bogs  along  the 
margins  of  the  former  glacial  lake.  Peat  accumulation  continued  until 
approximately  1500  years  B.P.  Wood  remains  in  the  Sphagnum  peat  are 
very  scarse,  indicating  a  bog  with  only  scattered  trees. 

The  change  from  Sphagnum  peat  to  a  nanolignid  forest  peat  marks 
the  raising  of  the  bog  surface  with  the  aggradation  of  permafrost. 

The  lack  of  permafrost  development  in  this  area  prior  to  approximately 
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4000  years  B.P.  (Zoltai  and  Tarnocai,  1973;  Zoltai  and  Pettapiece, 
1978),  in  spite  of  sufficiently  cool  mean  annual  temperatures  (see 
Chapter  3)»  was  because  of  the  lack  of  suitable  sites.  Permafrost 
will  not  develop  in  ponded  fens  due  to  the  vast  reservoir  of  latent 
heat  and  high  thermal  transfer  rates.  By  contrast,  Sphagnum  sites  are 
very  conducive  to  permafrost  aggradation  because  of  their  unique  ther¬ 
mal  qualities  (Brown,  1966) .  The  Picea-Ledum-Rubus-Cladonia  community 
that  followed  permafrost  development  has  remained  to  the  present  day. 

The  unfrozen  Sphagnum  bogs  contained  within  the  paisa  complex  are 
the  result  of  permafrost  degradation.  Stratigraphically,  two  Sphagnum 
layers  are  divided  by  a  thin  horizon  of  forest  peat  (Figure  4);  radio¬ 
carbon  dating  of  the  forest  peat  from  core  E  11,  at  a  depth  of  1.5  m, 
indicates  initial  freezing  took  place  approximately  3500  years  B.P. 

(WAT  400) .  The  areas  around  the  fringes  of  the  bogs  have  recently  re¬ 
frozen,  permafrost  being  re-established  approximately  400  years  B.P. 

The  marginal  areas  of  the  thaw  bogs  appear  to  be  comparatively  dynamic; 
continued  climatic  warming  may  result  in  another  collapse,  alterna¬ 
tively  cooling  may  re-establish  permafrost  along  the  bog  edges. 

The  *  thermokarst  *  sinkholes  by  contrast,  do  not  appear  to  be  a  re¬ 
sult  of  permafrost  degradation.  Sinkholes  differ  from  the  unfrozen 
bogs  in  that  they  are  contained  ponds  with  virtually  no  Sphagnum  and 
very  little  aquatic  vegetation  ( see  Chapter  1  and  Figure  4) .  The  peat 
stratigraphy  of  the  sinkholes  is  greatly  different  from  the  rest  of  the 
paisa  complex,  consisting  of  25  cm  of  moderately  decomposed  duckweed, 

15  cm  of  gyttja,  25  cm  of  ponded  fine  grained  detrital  sands,  and 
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glaciolacustrine  clay.  There  is  no  evidence  of  collapsed  peat  in  the 
center  of  the  hole.  Around  the  edges  of  the  sinkhole  however,  large 
lobes  of  semi -submerged  Sphagnum  peat  are  evident,  as  are  actively 
collapsing  walls.  It  is  unlikely  that  decomposition  rates  have  been 
sufficient  to  totally  digest  the  collapsed  peat.  The  stratigraphy 
suggests  that  some  holes  are  primary  features  and  have  been  present 
since  the  initial  aggradation  of  ground  ice.  Recent  permafrost  degra¬ 
dation  in  the  form  of  wall  collapse  has  enlarged  the  area  of  the  ponds. 

In  summary,  the  following  sequence  of  landscapes  have  developed 
in  peatlands  near  Ft.  Simpson: 

1.  a  large  deep  glacial  lake 

2.  a  partially  dessicated  lake  plain  with  numerous  shallow  ponds 

3.  a  paludificated  fenland 

4.  ombrotrophic  bogs,  and 

5.  extensive  permafrozen  peatlands. 


CHAPTER  5:  PERMAFROST  AGGRADATION  AND  DEGRADATION 
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5.1  Introduction 

The  development  of  ground  ice  in  peatlands  is  a  complex  process. 
While  cold  air  temperatures  are  a  necessary  criteria  for  permafrost 
aggradation,  paleotemperature  analysis  has  shown  that  it  is  not  a 
solely  sufficient  control.  The  role  of  different  vegetation  commun¬ 
ities  in  creating  various  soil  temperatures  has  "been  emphasized  by 
numerous  researchers  (Williams,  1968;  Brown,  1970).  Similarly,  the 
role  of  surface  water  in  modifying  ground  thermal  regimes  is  well 
acknowledged  (Mackay,  1972;  Brown,  197 0).  A  less  discussed  aspect 
of  the  role  of  vegetation  and  water,  however,  is  the  evolution  of  the 
vegetative  communities  and  the  hydrological  regimes  with  their  sub¬ 
sequent  effect  on  ground  ice  development.  This  evolution  is  dis¬ 
cussed  in  this  chapter.  An  examination  will  also  be  made  of  the  role 
of  the  mineral  soil  morphology  and  stratigraphy  in  determining  where 
in  a  watershed  permafrost  will  develop. 

The  melting  or  degradation  of  permafrost  is  equally  complex  as 
its  aggradation.  Three  components  of  the  melt  process:  lateral  melt¬ 
ing,  surface  warming,  and  melting  from  the  bottom  contribute  to  ground 
ice  degradation.  The  relative  roles  of  each  of  these  components  is 
examined  in  the  second  section  of  this  chapter.  The  rate  of  lateral 
melting  is  determined  quantitatively  from  historical  records.  By  ex¬ 
trapolation  of  the  rate  of  lateral  degradation,  the  former  extent  of 
permafrost  in  the  peatland  is  estimated. 


Both  ground  ice  aggradation  and  degradation  have  created  the 


present  day  paisa  relief.  This  fact  has  not  been  acknowledged  in  the 
literature.  Only  the  role  of  segregated  ground  ice  development  in  the 
mineral  soil  has  been  cited  as  a  contributor  to  paisa  height.  The  last 
section  of  this  chapter  shows  that  degradational  processes  are  equally 
important  in  creating  the  present  day  landscape. 

5.2  Ground  Ice  Aggradation 

5.21  Processes  and  Controls 

5.211  Climatic  Controls 

Two  climatic  requirements  for  the  presence  of  paisas  are 
a  Mean  Annual  Temperature  colder  than  -1°  C.  to  -3°  0.  (Ruuhijarvi, 
i960)  and  a  winter  precipitation  less  than  30  cm  (Lundqvist,  1962). 

The  present  day  climate  at  Fort  Simpson  is  within  these  bounds,  with  a 
Mean  Annual  Temperature  of  -4°  C.  and  a  winter  precipitation  averaging 
15  cm.  Holocene  Mean  Annual  Temperatures  at  Fort  Simpson  were  also  be¬ 
low  these  limits,  varying  between  -7°  C.  and  -3°  C.  (see  Chapter  3) • 

While  it  is  problematic  whether  paisas  are  only  being  maintained 
at  the  above  temperatures  or  whether  new  permafrost  is  actually  aggrad¬ 
ing,  it  is  still  unlikely  that  climate  has  been  an  impediment  to  perma¬ 
frost  development  over  the  past  8000  years.  The  question  remains  there¬ 
fore,  of  why  extensive  permafrost  did  not  develop  until  approximately 
4000  years  B.P.  in  spite  of  a  sufficiently  cold  climate?  Clearly  ed- 


aphic  factors  such  as  peatland  vegetation  and  hydrology  have  played  a 
key  role. 


5.212  Vegetational  And  Hydrological  Controls 

The  influence  of  vegetation,  in  particular  Sphagnum  moss, 
in  creating  ground  thermal  conditions  suitable  for  permafrost  aggrada¬ 
tion  has  "been  well  documented  (Brown,  1963;  Williams,  19 68;  Zoltai 
and  Tarnocai,  1975).  The  thermal  conductivity  of  Sphagnum  changes  suf¬ 
ficiently  with  moisture  content  or  when  frozen  that  the  outward  trans¬ 
fer  of  heat  in  the  winter  exceeds  the  inward  heat  transfer  in  the  sum¬ 
mer,  hence  creating  a  negative  heat  budget.  Because  of  this  unique 
thermal  characteristic,  permafrost  will  develop  beneath  Sphagnum  at  a 
higher  Mean  Annual  Temperature  than  beneath  mineral  soil. 

Nearly  all  paisas  have  a  Sphagnum  cover  at  the  time  of  formation; 
this  is  not  the  only  prerequisite  for  permafrost  aggradation  however. 
Paisas  from  the  study  area  typically  have  a  .8  meter  cap  of  Sphagnum, 
representing  approximately  1000  years  of  accumulation,  before  perma¬ 
frost  developed.  Also,  modern  collapse  bogs  contained  within  the  paisa 
complex  have  a  Sphagnum  cover  up  to  1.5  meters  in  thickness,  yet  are 
unfrozen.  This  dilemma  may  be  due  to  the  changing  hydrology  of  the 
bogs.  If  the  bogs  have  very  high  water  tables,  the  moss  will  be  con¬ 
tinually  saturated.  The  thermal  properties  of  the  Sphagnum  will  not 
seasonally  change;  a  negative  heat  budget  will  not  develop.  A  domed 
or  raised  bog  surface,  that  allows  partial  drying  of  the  upper  peat 
layers,  is  essential  to  permafrost  development  in  the  present  climate. 
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Black  Spruce  may  also  play  an  important  role  in  permafrost  aggra¬ 
dation  in  peatlands.  Small  permafrost  lenses  have  been  observed  be¬ 
neath  spruce  seedlings  that  were  growing  on  isolated  Sphagnum  hummocks 
in  an  otherwise  unfrozen  bog  (Zoltai  and  Tamocai,  1975)*  Shading  from 
the  spruce  creates  a  cooler  micro-climate  suitable  for  ground  ice  dev¬ 
elopment.  Zoltai  and  Tarnocai  maintain  that  the  small  frozen  hummocks 
eventually  coalesce  to  create  a  peat  plateau.  Further,  the  presence 
of  these  frozen  hummocks  indicates  that  paisas  are  presently  forming 
(Zoltai  and  Tarnocai,  1975). 

While  a  number  of  'embryonic'  paisas  of  the  type  described  were 
observed  near  the  study  site,  they  appear  to  be  temporary  features. 

One  frost  mound,  first  observed  in  August,  1976,  was  about  35  cm  high, 
had  a  basal  diameter  of  1.1m  with  ground  ice  to  a  depth  of  50  cm.  Two 
8  year  old  Black  Spruce  trees  grew  on  the  frost  mounds.  In  September 
of  the  following  year  the  mound  was  only  20  cm  high,  had  a  basal  diam¬ 
eter  of  0.6  m  with  ground  ice  at  a  depth  of  35  cm.  Two  other  frost 
mounds,  only  20  to  25  cm  high,  with  a  basal  diameter  of  30  cm,  in  the 
summer  of  1976,  had  completely  disappeared  by  September,  1977.  In  addi¬ 
tion,  a  ground  ice  lens  had  developed  beneath  a  nearby  seedling  covered 
hummock  where  the  previous  year  there  had  been  no  permafrost. 

Previous  discussion  on  the  observed  thermal  structure  of  the  paisa 
(Chapter  l)  has  demonstrated  the  sensitivity  of  the  permafrost  to  small 
changes  in  surface  temperatures.  Ruuhijarvi  (i960),  Wramner  (1967, 
1973)  and  Sollid  and  Sorbel  ( 197*0  have  also  noted  the  temporary  nature 
of  small  frost  mounds  in  peatlands. 
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The  strongest  argument  against  the  role  of  spruce  in  paisa  dev- 
eolopment  is  the  lack  of  spruce  on  paisas  in  Scandinavia  and  Ungava. 
Here,  paisas  develop  readily  above  the  tree  line. 

In  summary,  the  presence  of  large  ombro trophic  Sphagnum  bogs  that 
have  sufficiently  raised  and  drained  surfaces  are  necessary  for  ex¬ 
tensive  permafrost  development  in  peatlands  in  the  Ft.  Simpson  area. 

The  ombrotrophic  conditions  are  created  by  a  reorganization  of  the 
surface  drainage,  as  was  detailed  in  the  previous  chapter.  The  role 
of  Black  Spruce  is  important  in  maintaining  permafrost,  but  probably 
has  little  to  do  with  regional  permafrost  aggradation. 

5.213  Influence  of  Mineral  Soil  Topography 

The  possible  role  of  the  mineral  soil  topography  in  de¬ 
termining  preferential  zones  of  permafrost  aggradation  was  first  sug¬ 
gested  by  Thie  (197^) •  In  the  Ft.  Simpson  area,  the  close  association 
between  beach  ridges  and  paisa  complexes  (see  Chapter  l)  suggests  that 
the  initial  mineral  soil  topography  may  affect  the  sites  of  segregated 
ground  ice  development  in  the  peatlands.  The  topography  of  the  mineral 
soil  beneath  a  paisa  can  be  re-created  by  thawing  a  paisa  core  in  in¬ 
crements  under  in  situ  overburden  stresses.  Any  ’residual  height'  fol¬ 
lowing  thaw  consolidation  will  theoretically  equal  the  relief  of  the 
mineral  soil  prior  to  freezing. 

Methods  for  such  tests  have  recently  been  proven  for  a  variety  of 
soil  types  (Watson  et  al,  1973;  Luscher  and  Afifi,  1973;  Smith  and 
Morgenstern,  1972)  and  a  comprehensive  theory  for  thaw  consolidation 
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established  (Morgenstern  and  Nixon,  1971;  Nixon  and  Morgenstern,  1973)* 


The  total  thaw  consolidation,  E°,  of  a  thawing  soil  is  made  up  of 
three  separate  processes: 

1.  settlement  due  to  volume  reduction  associated  with  a  phase  change 
of  ice  to  water 

2.  consolidation  during  the  period  of  thaw 

3.  consolidation  subsequent  to  thaw 
The  total  thaw  strain  can  be  expressed: 

E°  =  V  +  ao  Jo  (p  +  x)  H 

where  E°  =  total  settlement 

Aq  =  initial  thaw  settlement  parameter 
aQ  =  coefficient  of  compressibility 
P  =  surcharge  load 
X  =  depth  to  the  thaw  front 
^ '  =  submerged  unit  weight  of  thawed  soil 

(Lusher  and  Afifi,  1973) 

Total  thaw  settlement  can  also  be  correlated  with  some  of  the  bulk 
properties  of  the  frozen  soil.  Speer  et  al  (1973)»  i-n  a  study  of 
ground  ice  variability  near  Norman  Wells,  found  an  excellent  correla¬ 
tion  between  frozen  bulk  density  and  total  thaw  strain  at  a  stress  eq¬ 
ual  to  the  average  stress  encountered  in  situ.  Following  this  method¬ 
ology,  31  thaw  consolidation  tests  were  made  on  the  frozen  mineral  soil 
from  core  E  3  of  the  paisa.  The  results  of  the  testing  are  tabulated 
in  Appendix  3*  The  testing  methods  were  similar  to  those  described  by 
Roggensack  (1978);  a  brief  summary  of  which  can  be  found  in  Appendix  4. 
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Following  consolidation,  linear  and  logarithmic  relationships  were 
derived  between  settlement  parameters  such  as  E°,  aQ,  and  Aq;  and  core 
physical  parameters  such  as  water  content  (w),  volumetric  ice  content 
(v),  and  frozen  hulk  density  (Figures  22  to  26  ).  The  most  signif¬ 

icant  of  these  correlations  are  listed  in  Table  8  .  The  E*  f  re- 
lationship  and  E°:  Aq  +  aQ  relationship  gave  the  most  accurate  results; 
both  predicted  settlements  within  1.5  %  of  that  measured. 

The  bulk  density  relationship  was  used  to  predict  settlement  in 
six  other  cores  from  the  paisa  (Appendix  3).  Total  thaw  settlements 
for  each  core  are  listed  in  Table  9  ,  as  are  the  respective  differences 
in  elevation  between  the  peat -mineral  soil  interface  in  the  paisa  and 
that  in  the  fen. 

In  each  of  the  cores,  a  positive  residual  mineral  soil  relief  was 
detected.  The  magnitude  of  the  relief,  from  .23  to  .65  m  represents 
from  9 %  to  4 6%  of  the  total  frozen  mineral  soil  relief,  well  in  excess 
of  the  methodological  error.  The  residual  relief  is  approximately  eq¬ 
ual  to  the  thickness  of  the  sand  layer  usually  encountered  immediately 
beneath  the  peat.  Freezing  of  the  sand  layer  will  produce  insignificant 
heave  compared  to  freezing  of  a  silty  clay  layer.  If  the  initial  soil 
relief  beneath  the  paisa  is  largely  provided  by  a  sand  hummock,  then 
the  total  heave  in  the  mineral  soil  will  decrease  with  the  relative 
thickness  of  sand  layer  that  is  present.  A  plot  of  total  heave  in  the 
mineral  soil  vs.  initial  mineral  soil  relief  shows  this  to  be  the  case 

(Figure  27),  with  a  very  strong  negative  correlation  between  the  two 
parameters . 
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THAW  SETTLEMENT  PARAMETER  (Ao) 


FIGURE  22 


-  A0  =  -6.32  x  I04  Jff  +  I  17.2 

- Ao  =  -658.6  -  104.74  In  / f 


THAW  SETTLEMENT  PARAMETER,  A0  ,  FOR  MINERAL 
SOIL  RELATED  TO  FROZEN  BULK  DENSITY. 
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TABLE  9 


TOTAL  THAW  SETTLEMENT  OF  THE  MINERAL  SOIL  AND 
OBSERVED  FROZEN  MINERAL  SOIL  RELIEF  IN  THE  PALSA  COMPLEX 


Residual  mineral 


Core 

Thaw  settlement 
(m) 

Observed  frozen  soil 
relief  (m) 

soil  relief 
thaw  (m) 

S  100  * 

1.55 

1.95 

+  .40 

S  101 

1.38 

1.93 

+  .35 

S  78 

1.0  5 

1.52 

+  .47 

S  110 

0.93 

1.49 

+  .58 

E  3 

2.54 

2.76 

+  .23 

E  11 

0.77 

1.42 

+  .65 

Thaw  settlements  "based  on:  E°  =  -  .  IO36  ^  +  182.7 


* 


Thaw  settlement  directly  measured  in  Permode 


. 
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Further  evidence  of  the  hummocky  mineral  soil  interface  that  under¬ 
lies  the  paisa  complex  is  the  irregular  mineral  soil  elevations  meas¬ 
ured  in  the  collapse  hogs  and  sinkholes  (Table  10).  Differences  of  up 
to  4  m  exist  between  sites,  but  the  average  differential  relief  is  ap¬ 
proximately  .5  to  1.0  meters. 

Changes  in  the  soil  consolidation  properties  due  to  freezing  will 
not  significantly  affect  the  above  results.  Local  overconsolidation 
of  the  mineral  soil  can  occur  during  freeze-thaw  cycles  due  to  negative 
pore  pressures  generated  during  freezing  (Smith,  1972;  Nixon  and  Morg- 
enstern,  1973)*  Indeed,  residual  stress  determinations  (Table  2,  Ap¬ 
pendix  3 )  on  the  thawed  permafrost  core  show  significant  residual 
stress  in  the  massive  silts,  indicating  local  overconsolidation.  Be¬ 
cause  of  soil  discontinuities  due  to  lensing  however,  the  bulk  soil 
porosity  of  a  soil  after  freezing  and  thawing  often  increases  over  the 
original  (Tsytovich,  1976)  especially  under  low  effective  stresses. 

These  two  processes,  local  ped  overconsolidation  and  an  overall  bulk 
dilatancy  due  to  lensing,  work  in  opposite  directions  in  determining 
the  final  thaw  settlement;  but  the  most  characteristic  phenomenon  will 
be  an  increase  in  the  bulk  void  ratio  of  the  thawed  soil  (Tsytovich, 
1976).  Rates  of  settlement  and  pore  pressures  in  the  thawing  soil 
will  consequently  be  affected,  however  changes  in  the  total  thaw  set¬ 
tlement  due  to  these  processes  will  be  minimal. 

In  conclusion,  the  paisa  complex  has  preferentially  developed  over 
an  area  with  a  slightly  raised,  irregular,  mineral  soil  surface.  This 
mineral  soil  relief  contributes  from  10  to  25  per  cent  of  the  observed 


height. 
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5.3  Ground.  Ice  Degradation 


Restricted  to  the  discontinuous  permafrost  zone,  a  paisa  is  in 
a  delicate  thermal  balance  with  its  physical  environment.  A  shift  in 
this  balance  can  produce  ground  ice  degradation.  Many  degradational 
processes  can  be  attributed  to  climatic  change,  but  others  are  caused 
by  changes  in  edaphic  factors.  This  section  examines  ground  ice  deg¬ 
radational  processes  operating  on  the  paisa  complex,  determines  rates 
of  permafrost  degradation  and  attempts  to  determine  the  former  extent 
of  permafrost  in  the  area. 

5.31  Processes  and  Controls 

5.311  Surface  Processes 

Surface  cracking  has  been  considered  by  some  authors  to 
be  the  most  important  degradational  process  (Ruuhijarvi,  i960;  Sven- 
sson,  1962) .  The  surface  peat  layer  of  the  study  paisa  is  very  heav¬ 
ily  cracked.  Crack  dimesions  run  up  to  6  m  in  length,  varying  from  a 
few  millimeters  to  80  cm  in  width.  Cracking  is  most  extreme  on  the 
paisa  summits,  on  the  steeper  slopes,  and  particularly  near  the  col¬ 
lapsing  edges  where  there  is  a  strong  crack  orientation  parallel  to 
the  edges . 

The  cracks  are  the  result  of  a  number  of  possible  causes: 

1.  surface  rupture  due  to  differential  heaving  during  ground  ice  ag¬ 
gradation 


Partial  infilling  of  the  cracks  with  semi -decomposed,  vegetation,  in 
particular  Ledum  leaves,  may  partially  offset  their  thermal  influence. 
Only  in  heavy  rainstorms  does  convective  heat  increase  the  depth  of 
the  crack  relative  to  the  active  layer.  This  increase  was  always 
observed  to  disappear  within  a  few  days  of  cessation  of  the  rain. 

There  is  no  evidence  to  suggest  that  the  present  day  active  layer  is 
increasing . 

5. 312  Bottom  Melting 

Bottom  melting  of  a  paisa  has  been  discussed  in  detail 
in  Chapter  One.  It  was  shown  that  a  mean  annual  surface  temperature 
at  the  study  site  greater  than  -0.75°  0.  would  result  in  melting  at 
the  base  of  the  permafrost  at  a  rate  of  approximately  1  cm/year.  This 
melting  rate  is  entirely  governed  by  the  geothermal  heat  flux,  but  can 
vary  depending  on  the  local  ice  stratigraphy. 

5.313  Lateral  Erosion 

Lateral  erosion  of  the  paisa  complex  occurs  mainly 
through  thermal  degradation  of  the  paisa  perimeter  by  ponded  water. 

The  process  is  facilitated  by  the  high  thermal  conductivity  of  sat¬ 
urated  peat  as  compared  to  the  dry  peat  on  the  paisa  summit  as  well 
as  by  the  high  convective  heat  capacity  of  water. 


Where  the  paisa  sides  are  in  direct  contact  with  water,  the 
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permafrost  core  is  undermined,  creating  a  steep  erosional  edge. 
Undercutting  of  the  peat  hank  is  greatest  at  a  depth  of  approximately 
30  to  50  cm  beneath  the  pond  surface,  corresponding  to  the  thickness 
of  the  winter  ice  layer.  The  ice  seasonally  protects  part  of  the 
paisa  sides;  below  this  layer  the  submerged  permafrost  is  being 
continuously  eroded.  The  result  is  a  permafrost  'nose*  near  and 
just  above  the  water  table  (Figure  4  ).  Zoltai  and  Tarnocai  (1971) 
have  also  noted  such  features  but  have  considered  them  proof  of  ground 
ice  aggradation  instead  of  degradation. 

Numerous  surface  cracks  form  parallel  to  the  collapsing  edge  of 
the  paisa.  The  cracking  results  in  discrete  peat  blocks  that  slide 
along  the  permafrost  surface  when  sufficient  undercutting  has  occurred. 
Progressive,  widening  of  the  cracks  toward  the  collapsing  edge  in¬ 
dicates  an  accelerating  creep  rate  towards  the  margin  (Plate  2  ). 

Lateral  erosion  may  be  accompanied  by  concurrent  frost  heave 
along  the  paisa  perimeter.  Most  paisa  complexes  and  peat  plateaus 
have  distinct  but  discontinuous  marginal  ridges,  usually  from  .5  m  to 
1.0  m  above  the  interior  regions  of  the  plateau.  As  ponded  water  is 
available  along  most  of  the  collapsing  perimeter,  water  may  be  mig¬ 
rating  into  the  frozen  soils  along  thermal  gradients,  producing 
'secondary  heave'  in  the  paisa  soils.  Further  evidence  for,  and  dis¬ 
cussion  of,  secondary  heave  is  in  Chapter  Six. 


The  paisa  complex  also  erodes  laterally  from  interior  sinkholes 
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2.  thermal  contraction  of  the  surface  layers,  and 

3.  dessication  cracking. 

Cracks  up  to  several  meters  in  length,  hut  only  a  few  mil¬ 
limeters  wide  were  observed,  suggesting  that  cracking  starts  while 
the  peat  is  frozen.  Clearing  of  the  snow  from  a  paisa  summit  dur¬ 
ing  the  winter  revealed  that  some  of  the  thin  cracks  were  infilled 
with  ice.  Thermal  gradients  due  to  cold  penetration  in  the  cracks 
may  induce  moisture  migration  towards  the  crack  from  the  adjacent 
peat.  Crack  enlargement  due  to  subsequent  ice  wedging  is  con¬ 
sidered  a  minor  process  however,  and  limited  to  the  narrowest  cracks. 

Dessication  cracking  may  be  a  relatively  more  important  process 
on  the  steeper  sloping  paisa  sides.  Rainfall  and  snowmelt  quickly 
run  off  the  steep  slopes  with  little  water  infiltration,  producing 
a  locally  dessicated  environment.  The  only  vegetation  growing  on  the 
slopes  is  Ledum,  which  grows  within  the  wider  cracks,  presumably  due 
to  the  higher  moisture  availability. 

Overall,  surface  melting  is  not  considered  to  be  a  major  deg- 
radational  process  at  the  study  site.  Dry  Sphagnum  peat  has  a 
thermal  conductivity  as  low  as  .05  Cal/m-hr.-°C  (Jumikis,  1977), 
providing  excellent  insulation  from  surface  warming.  Surface  cracks, 
while  extensive,  do  not  have  a  strong  effect  on  the  overall  surface 
insulation.  The  average  depth  to  permafrost  beneath  the  cracks  does 
not  exceed  the  average  active  layer  depth  on  other  parts  of  the  paisa. 


. 


e:  X' -O  git'  -i  :  rers-T  ij.-C|L  <_  v-  i  *•£  qj.  istr.  it  -do 


■ 


and  collapse  bogs.  The  mechanism  of  erosion  is  similar  to  erosion 
along  the  outside  perimeter,  but  involves  much  slower  degradation  of 
the  ground  ice  without  a  significant  collapsing  edge.  The  bogs  are 
filled  with  Sphagnum,  greatly  reducing  rates  of  water  filtration  and 
hence  heat  transfer.  Further,  the  near  surface  summer  water  temp¬ 
eratures  in  the  bog  areas  average  only  1.5°  C.  compared  to  5  to  7°  C. 
in  the  surrounding  fens;  the  lower  temperatures  being  due  to  a  comb¬ 
ination  of  the  high  evapotranspiration,  high  albedo,  and  low  thermal 
conductivity  of  the  Sphagnum  moss. 

Of  the  three  degradational  processes,  lateral  erosion  is  by  far 
the  most  effective.  This  is  borne  out  by  the  extraordinary  rates  of 
lateral  melting  that  have  occurred  over  the  past  30  years,  as  is  dis¬ 
cussed  in  the  following  sections. 

5.32  Former  Extent  of  Permafrost 

The  former  extent  of  permafrost  near  the  study  site  can  be  in¬ 
ferred  from  the  present  day  distribution  of  peatland  types.  As  pre¬ 
viously  discussed  (Chaper  l),  the  peatland  type  described  as  "open 
sedge  fen  with  abundant  non-oriented  ponds"  preferentially  occurs 
in  areas  known  to  have  been  formerly  permafrozen.  Cores  taken  from 
this  peatland  type  should  therefore  stratigraphically  indicate  the 
former  presence  of  a  thick  Sphagnum  layer  and  possibly  a  Black  Spruce 
forest  cover.  These  characteristics  were  indeed  found  in  a  core  ob¬ 
tained  from  a  sedge  fen  1.2  km  north  of  the  paisa  complex. 
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Chemical  analysis  was  used  to  detect  the  Sphagnum  because  the 
plants  were  too  decomposed  to  identify  macroscopically .  Major  cation 
differences  occur  in  the  vegetation  from  each  peatland  type  because 
of  their  distinctive  aquatic  nutrient  regimes  (Stanek  et  al,  1977). 
While  decomposition  processes  will  modify  the  peat  chemistry,  changing 
particularly  the  carbon  and  nitrogen  contents,  cation  differences  are 
still  evident,  hence  providing  a  stratigraphic  tool  for  peat  differen¬ 
tiation. 

Results  from  the  chemical  analysis  are  presented  in  Table  11. 
Samples  from  modern  peat-forming  plants  indicate  that  Ca  is  the  most 
sensitive  indicator  for  distinguishing  peat  types.  Sphagnum  mosses 
and  forest  peat  typically  have  Ca  contents  of  2  to  6  mg/gm  while  fens 
average  10  to  22  mg/gm.  The  chemical  differences  are  maintained  in 
spite  of  moderate  decomposition;  proof  of  which  is  seen  in  the  re¬ 
sults  from  the  analysis  of  the  frozen  paisa  core.  The  peat  in  this 
core  is  identifiable  macroscopically;  the  stratigraphy  of  which  cor¬ 
relates  well  with  Ca  differences  between  the  ombrotrophic  and  minera- 
trophic  peats  (Table  11 ) .  The  peat  core  from  the  fenland  shows  a 
shift  from  a  high  Ca  surface  peat  to  a  low  Ca  peat  at  a  depth  of  40 
to  70  cm,  then  a  very  high  Ca  lower  peat.  This  is  interpreted  as  a 
basal  fen  peat,  a  middle  Sphagnum  peat  followed  by  a  return  to  minera- 
trophic  fen  conditions.  A  significant  period  of  ombrotrophic  con¬ 
ditions  once  existed  at  this  site. 


The  presence  of  Black  Spruce  fragments  in  the  peat  also  suggests 
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a  previously  frozen  peatland,  as  modern  spruce  are  nearly  always  con¬ 
fined  to  permafrozen  areas  in  the  peatland.  Abundant  pieces  of  wood 
charcoal,  less  than  2  mm  in  size,  were  observed  along  a  horizon  at  a 
depth  of  43  cm  in  the  peat  core.  This  horizon  is  near  the  top  of  the 
Sphagnum  zone  as  defined  by  chemical  analysis,  suggesting  that  a 
former  Black  Spruce  forest  may  have  been  the  source.  D.  Rennie 
(personal  communication)  has  also  found  identifiable  pieces  of  Black 
Spruce  in  a  number  of  cores  taken  from  the  fen  in  this  vicinity. 

The  former  presence  of  Sphagnum  and  Black  Spruce  communities  in 
the  "disoriented  pond-sedge  fen"  peatland  supports  the  hypothesis  that 
this  peatland  type  is  an  indicator  of  previously  permafrozen  terrain. 
As  such,  the  study  site  was  once  part  of  a  large  permafrost  body  cov¬ 
ering  approximately  600  hectares.  The  present  day  paisa  and  peat 
plateau  are  erosional  remnants  of  the  once  more  extensive  landform. 
Permafrost  is  not  thought,  however,  to  have  been  continuous.  Most  of 
the  mineratrophic  fens  near  the  large  lakes  in  the  middle  of  the  lake 
basin,  to  the  north  of  the  study  site,  do  not  display  the  characteris¬ 
tic  peatland  type  of  previously  frozen  areas. 
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5* 33  Hates  of  Ground  Ice  Degradation 

Rates  of  permafrost  degradation  can  be  measured  by  comparing 
recent  and  older  aerial  photography.  This  approach  has  been  used  by 
Thie  (197*0  in  northern  Manitoba,  and  by  Reid  (1977)  in  northern 
Alberta.  Aerial  photography  for  the  study-  site  was  available  from  1948 
(1:40,000  RCAF)  and  1976  (1:25,000,  Burnett  Resources).  Both  the  lin¬ 
ear  rate  of  retreat  of  the  paisa  complex  margin  and  the  rate  of  areal 
degradation  were  measured. 

To  measure  changes  in  the  areal  extent  of  the  permafrost,  indiv¬ 
idual  peat  plateaus  were  traced  onto  'Mylar'  film,  carefully  cut  out 

-4 

with  a  dissecting  scalpel,  then  weighed  to  an  accuracy  of  1  x  10  gm. 
Mylar  film  has  extremely  consistent  density  and  thickness,  permitting 
conversion  from  weight  to  area.  Paisas  and  peat  plateaus  greater  than 
200  m  in  diameter  have  degraded  by  15  to  20 %  over  the  last  28  years. 
Medium  sized  plateaus  with  a  1948  diameter  of  50  to  150  meters  exper¬ 
ienced  considerably  faster  decay,  losing  an  average  of  44#  of  their 
1948  area.  Most  permafrost  islands  having  a  diameter  of  less  than  50  m 

in.  1948  were  no  longer  discernible  in  1976  photography. 

/ 

Rates  of  margin  retreat  were  measured  relative  to  a  line  perp¬ 
endicular  to  the  1976  paisa  edge.  A  magnifying  stereoscope,  contain¬ 
ing  a  built  in  scale  with  gradations  of  0.1  mm,  was  used  for  taking 
measurements.  The  rate  of  edge  retreat  varied  considerably  around 
the  paisa:  from  8  to  38  meters  for  paisas  in  excess  of  200  meters; 
from  21  to  43  meters  for  paisas  of  50  to  150  meter  diameter.  No 


. 


' 
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particular  aspect  showed  preferential  melting.  The  size  of  thermo¬ 
karst  holes  and  collapse  hogs  in  the  middle  of  the  paisa  complex  was 
relatively  unchanged  from  1948  to  1976.  These  rates  of  edge  col¬ 
lapse  are  approximately  50 %  less  than  that  measured  by  Reid  (1977) 
near  Zama  Lake,  Alberta  and  are  approximately  the  same  as  that  re¬ 
ported  by  Thie  (1974)  in  northern  Manitoba. 


On  the  basis  of  the  former  extent  of  the  permafrost,  it  appears 
that  rates  of  degradation  have  been  greater  than  aggradation  for 
a  considerable  number  of  years.  If  the  former  size  of  the  paisa 
complex  has  degraded  at  present  day  melting  rates,  then  all  the 
degradation  may  have  occurred  in  the  past  250  years. 


5.4  The  Effect  on  Paisa  Relief 

Paisas  have  been  traditionally  viewed  as  ground  ice  aggradational 
features;  their  height  measured,  using  the  surrounding  fen  as  the 
local  base  level.  The  relief  is  considered  a  result  of  segregational 
ice  development  in  the  mineral  soil.  As  segregated  ice  lenses  are 
rarely  observed  in  frozen  peat,  the  peat  thickness  in  the  unfrozen  fen 


■ 
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and  on  the  frozen  paisa  are  assumed  identical,  hence  contribute  little 
to  the  paisa  relief. 

This  section  presents  evidence  that  contradicts  this  aggrada- 
tional  theory  of  paisa  development.  The  relative  roles  that  degrada- 
tional  and  aggradational  processes  play  in- the  resultant  relief  of  the 
paisa  complex  is  examined. 

Various  lines  of  evidence  have  shown  that  the  area  surrounding 
the  study  site  has  recently  thawed.  A  nearly  continuous  hand  of 
frozen  peatland  may  have  once  existed  along  the  margins  of  the  former 
glacial  lake,  as  late  as  400  years  Before  Present.  With  recent  clim¬ 
atic  warming,  ground  ice  degradation  has  occurred  at  an  extremely  rap¬ 
id  rate;  the  modern  paisa  complex  is  a  remnant  of  the  former  perma¬ 
frost  landform.  The  surrounding  unfrozen  peatland  is  composed  of  both 
the  thaw  consolidated  peat  and  the  thaw  consolidated  mineral  soil. 

Thaw  consolidation  tests,  as  detailed  above,  recorded  an  average 
settlement  of  1.3  meters  in  the  mineral  soil  of  the  paisas  following 
thaw.  This  settlement  accounts  for  less  than  fifty  percent  of  the 
average  3.0  meter  paisa  height  above  the  fen.  As  there  is  no  evidence 
for  wildly  fluctuating  fen  levels,  this  conclusion  clearly  shows  that 
segregational  ice  in  the  mineral  soil  is  not  solely  responsible  for 
paisa  relief. 


To  determine  what  was  responsible  for  the  remaining  height,  thaw 


. 


consolidation  tests  were  made  on  the  frozen  peat  portion  of  cores 
S  100  and  E  3.  The  consolidation  methods  were  identical  to  that  used 
on  the  frozen  mineral  soil.  A  total  of  18  samples  were  tested;  the 
results  are  listed  in  Appendix  3 .  Empirical  relationships  between  the 
thaw  settlement  and  physical  parameters  such  as  moisture  content  and 
bulk  density  were  not  as  precise  as  those  for  the  mineral  soils 
(Figures  28  to  29);  however,  the  simple  relationship  of  total  strain 
equal  to  the  volumetric  ice  content  of  the  frozen  peat  had  a  signif¬ 
icantly  high  correlation  (Figure  29).  This  relationship  predicted  the 
total  settlement  of  the  thawed  peat  cores  to  within  6.0  %  of  that 
measured  (Table  12). 

Application  of  the  thaw  settlement  equation:  E  =  V  +  4.5»  to 
the  five  other  frozen  peat  cores  shows  that  the  total  thaw  settlement 
of  the  peat  is  slightly  greater  than  the  thaw  settlement  of  the  min¬ 
eral  soils  (Table  13).  Further,  it  is  interesting  to  note  that  the 
moisture  contents  of  the  frozen  peat  are  approximately  equal  to  the 
unfrozen  peat  in  nearby  bogs  (Appendix  3  )>  a  finding  consistant  with 
other  researchers  results  (Zoltai  and  Tarnocai,  1975;  Reid  et  al, 
1977).  As  there  is  no  excess  water,  the  main  component  of  the  thaw 
consolidation  equation  (see  Section  5.213)  is  that  due  to  the  normal 
consolidation  of  the  peat  as  it  settles  under  its  own  weight.  The 
large  amount  of  observed  settlement  is  due  to  the  fact  that  the  peat 
must  have  frozen  while  still  in  the  unconsolidated  state.  When  melt¬ 
ing  finally  occurs,  the  peat  has  been  uplifted  above  the  fen  surface, 
hence  is  free  to  drain  and  consolidate.  A  combination  of  the  thaw  of 
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TABLE  12 

EMPIRICAL  RELATIONSHIPS  PREDICTING  THAW  SETTLEMENT  IN  PEAT 


CORRELATION  (r) 

1.  E°  =  -  .0162  Y  +  56.93  .02 

2.  E°  =  97  -  7.97  In  y  f  .05 

3.  E°  =  1.09  wy  +  4.36  .81 

OR 

E°  =  V  +  4.36  .81 


V  =  volumetric  ice  content  (%) 


TABLE  13 


THAW  SETTLEMENT  OF  THE  FROZEN  PEAT  GORES, 
AND  COMPARISON  WITH  MINERAL  SOIL  SETTLEMENT 


Frozen  peat 
Gore  thickness  ( m) 


Peat  thaw  Mineral  soil  thaw 

settlement  (m)  settlement  (m) 


S  100 

3.35 

1.79 

1.55 

S  101 

3.37 

1.89 

1.58 

S  78 

3.26 

2.46 

1.05 

s  no 

3.10 

2.51 

.93 

E  3 

3.20 

1.71 

2.54 

E  11 

4.30 

2.74 

.77 

Thaw  settlements  of  the  peat  based  on:  E°  = 
Measured  in  Permode 


V  +  4.5 


• 

segregated  ice  lenses  in  the  mineral  soil  plus  the  normal  consolida¬ 
tion  of  the  peat,  is  approximately  equal  to  the  amount  that  the  land 
surface,  now  covered  hy  the  fen,  has  dropped  below  the  former  perma¬ 
frost  surface.  The  permafrost  that  remains  forms  the  paisa  complex. 


■ 


CHAPTER  6:  CONCLUDING  REMARKS 
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6.1  Regional  Permafrost  Development  in  the  Fort  Simpson  Peatlands 

Permafrost  distribution  in  the  Fort  Simpson  area  is  very  dis¬ 
continuous.  It  is  found  sporadically,  to  various  depths,  in  treed 
till  uplands;  however,  the  largest  areal  extent  is  within  the  large 
glacial  lacustrine  basins.  Any  conceptual,  model  of  permafrost  dev¬ 
elopment  in  this  area  must  include  the  following  facts: 

1.  The  most  extensive  zones  of  segregated  ground  ice  accumulation  in 
the  Fort  Simpson  area,  occur  parallel  to  the  margins  of  the  former 
glacial  lakes  (Chapter  2). 

2.  These  areas  are  underlain  by  a  hummocky,  stratified  mineral  soil, 
slightly  elevated  above  the  surrounding  lacustrine  plain  (Chapter 

5). 

3.  Permafrost  has  developed  relatively  recently,  approximately  2000 
years  Before  Present  (Chapter  4). 

4.  This  late  development  was  in  spite  of  relatively  cold  temperatures 
throughout  the  Holocene  (Chapter  3) • 

5.  There  is  evidence  for  only  one  period  of  regional  permafrost  form¬ 
ation;  it  does  not  appear  to  have  been  forming  over  the  past  few 
hundred  years  (Chapter  5) • 

6.  Permafrost  development  was  always  preceded  by  an  accumulation  of 
Sphagnum  moss;  however,  as  much  as  1000  years  of  Sphagnum  could 
accumulate  before  permafrost  formed  ( Chapter  4) . 

There  is  evidence  for  a  period  of  aridity  sometime  during  the 
Holocene  (Chapter  2). 


7. 


« 


While  somewhat  conjectural,  the  following  scenario  best  fits 
the  findings  of  this  thesis.  During  the  waning  stages  of  Glacial 
Lake  McConnell,  lake  levels  fluctuated  greatly,  creating  extensive 
hummocky  beach  deposits  and  interfingered  sandy  and  silty  pro- 
lacustrine  sediments  along  its  margins.  When  the  lake  finally 
drained,  an  extensive  fenland  developed  over  the  lacustrine  basin. 

With  the  progressive  development  of  an  integrated  drainage  system, 
fen  water  levels  slowly  dropped.  Local  water  table  divides  developed 
over  the  more  elevated  parts  of  the  lake  basin,  especially  near  the 
beach  ridges.  The  water  divides  soon  became  ombrotrophic,  leading 
to  the  subsequent  development  of  extensive  Sphagnum  bogs. 

Permafrost  does  not  develop  in  saturated  Sphagnum;  it  is  only 
dry  Sphagnum  that  has  the  unique  insulative  properties  that  maintain 
cold  soil  temperatures.  A  period  of  aridity  approximately  2000  years 
B.P.  may  have  caused  sufficient  drying  of  the  Sphagnum  to  initiate 
permafrost  growth;  alternatively  water  levels  could  have  continued 
to  drop,  eventually  exposing  the  pro -lacustrine  zone.  The  temperature 
does  not  appear  to  have  been  a  factor,  remaining  relatively  cold 
throughout  the  Holocene.  Eventually  extensive  permafrost  develop¬ 
ment  occurred  over  the  bogs.  Evidence  suggests  that  permafrost  cov¬ 
ered  an  area  of  at  least  600  hectares  around  the  study  site.  Perma¬ 
frost  was  not  continuous  however,  as  there  is  evidence  that  the  deep¬ 
er  parts  of  the  lacustrine  basin  have  never  been  permafrozen. 


While  permafrost  in  the  till  uplands  may  be  as  deep  as  in  the 
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paisa  complex,  by  comparison  there  is  very  little  segregated  ground 
ice.  The  unique  stratigraphy  of  the  pro -lacustrine  environment,  with 
interfingering  layers  of  sand  and  silt,  provides  an  adequate  supply  of 
moisture  to  the  freezing  front,  ensuring  the  growth  of  large  ice 
lenses.  Chemical  and  isotopic  analysis  of  the  ground  ice  has  shown 
that  mineral  soil  groundwater,  not  the  ponded  fens,  is  the  water 
source  for  ice  lens  development,  hence  demonstrating  the  importance  of 
groundwater  flow  in  the  sand  lenses. 

The  ground  ice  stratigraphy  has  also  been  influenced  by  the  evol¬ 
ution  of  the  peatland.  A  characteristic  feature  of  nearly  all  cores 
is  a  large  ice  lens,  up  to  1.0  m  thick  at  the  peat -mineral  soil  inter¬ 
face.  This  thick  lens  may  be  partially  the  result  of  high  concen¬ 
trations  of  CaSO^  at  the  top  of  the  mineral  soil,  which  was  probably 
deposited  during  a  prior  dessication  period  (Chapter  2).  It  is  well 
known  that  solutes  will  depress  the  freezing  point  of  a  liquid,  the 
magnitude  of  which  is  proportional  to  the  solute  concentration  (Hallet, 
1978).  Experiments  by  R.D.  Miller  (personal  communication,  1978)  have 
shown  that  while  CaSO^  concentrations  of  100%  by  dry  weight  are  very 
effective  in  reducing  frost  heave,  concentrations  of  by  dry  weight 
have  virtually  no  effect  on  water  intake  rate.  The  net  effect  of  the 
solute  then  is  to  stall  the  advance  of  the  freezing  front  by  depres¬ 
sing  the  freezing  point  but  without  reducing  the  suction  potential. 

High  suction  potentials  are  ensured  by  the  very  decomposed  gyttja, 
with  a  large  internal  surface  area,  found  immediately  above  the  miner¬ 
al  soil.  The  combined  effect  produced  the  very  thick  ice  lens  at  the 


interface. 


. 
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In  summary,  the  location  of  permafrost  development,  its  ice 
content,  and  ice  stratigraphy,  have  all  been  influenced  by  prior  geo¬ 
logic,  climatic,  and  hydrologic  events  within  the  glacial  lacustrine 
basins  of  the  Fort  Simpson  area.  The  pro -lacustrine  zone  is  part¬ 
icularly  suitable  for  segregated  ground  ice  development. 

6.2  The  Paisa  Complex:  A  Relic  Thermo -erosional  Landform 

Historical  air  photo  evidence  has  proven  that  rapid  areal  perma¬ 
frost  degradation  is  occurring  around  the  study  site.  Stratigraphic 
chemical  results  indicate  that  large  areas  of  the  fenland  presently 
surrounding  the  modern  paisa  complex  were  once  permafrozen.  The 
evidence  is  strong  therefore  that  the  present  paisa  complex  -  fenland 
environment  is  the  degrading  remnant  of  a  once  extensive  permafrost 
landscape. 

Previous  experimental  results  and  discussion  have  demonstrated 
that  segregated  ice  development  in  the  mineral  soil  does  not  account 
for  the  total  paisa  relief  (Chapter  5) •  Segregated  ice  development 
in  the  mineral  soil  is  therefore  not  the  sole  cause  of  paisa  develop¬ 
ment.  Rather,  the  thaw  consolidation  of  the  frozen  peat  accounts  for 
a  major  part  of  the  observed  height  of  a  frozen  paisa  complex  above 
the  surrounding  unfrozen  fen.  Significantly,  the  main  component  of 
the  thaw  consolidation  of  the  peat  is  that  due  to  normal  consolidation. 
The  peat  had  frozen  in  a  buoyant,  unconsolidated  state  but  has  thawed 
while  uplifted,  hence  consolidating  under  its  own  saturated  weight 
(Chapter  5) • 


I 


Numerous  cores  from  the  fen  reveal  an  average  peat  thickness  of 
1.5  m>  considerably  less  than  that  of  the  frozen  peat.  While  some 
of  this  difference  can  be  ascribed  to  decomposition,  consolidation 
is  indicated  by  the  ubiquitous  presence  of  small  ponds  and  lakes  over 
the  thawed  zone  (Plate  l).  Such  a  concentration  of  ponds  is  not 
found  elsewhere  in  the  peatlands.  Svensson  (1969)  has  also  remarked 
that  circular  lakes  often  mark  sites  of  former  paisas  in  Northern 
Norway.  The  observation  of  a  water  moat  around  thawing  paisas  has 
also  been  made  in  numerous  case  studies  (Lundqvist,  1959;  Salmi, 
1970;  Svensson,  1969) .  With  aggradational  theories  of  paisa  devel¬ 
opment,  the  water  moats  have  been  ascribed  to  depressions  of  the  min¬ 
eral  soil  by  the  increasing  weight  of  the  growing  paisa  (Svensson, 
1969).  By  perceiving  paisas  as  degrading  landforms,  the  moats  are 
recognized  as  features  created  by  consolidation  of  the  thawing  peat. 

A  persistent  theory  in  paisa  development  is  the  notion  of  cyc¬ 
lical  aggradation  and  degradation  of  paisas.  There  is  no  evidence, 
in  the  Fort  Simpson  area,  for  more  than  the  present  occurrence  of 
paisas.  As  discussed  above,  ground  ice  development  is  the  net  result 
of  the  unique  evolution  of  the  peatland.  Permafrost  appears  to  have 
developed  during  a  single  interval  approximately  2000  -  3000  years 
B.P.,  when  there  was  the  right  combination  of  climate,  hydrology  and 
vegetation.  Over  the  past  few  hundred  years  it  has  been  continually 
degrading.  The  present  occurrence  of  the  paisa  complex  is  due  only 
to  the  insulative  properties  of  the  peat  cover.  Recognition  of  this 
will  necessitate  a  change  in  the  classification  of  paisas.  The  terms 
’young',  'mature',  and  'overmature'  (Zoltai  and  Tarnocai,  1975)  are 


« 


/» 


not  applicable  to  this  situation. 


It  is  recognized  that  some  paisas  are  in  fact  probably  aggrad- 
ational.  This  would  certainly  appear  to  be  the  case  for  small 
frozen  mounds  found  in  the  middle  of  Sphagnum  bogs  or  for  the 
'esker  paisas'  (Jahn,  1975)  found  along  streams  in  Northern  Europe. 
Certainly,  two  types  of  classification  are  needed.  On  an  areal 
basis  however,  the  degradational  type  of  paisa  is,  in  the  author's 
opinion,  the  most  prevalent. 

In  summary,  paisa  complexes  in  the  study  area  are  erosional 
remnants  of  a  once  much  larger,  thermally  degraded,  permafrost  land- 
form.  The  observed  relief  is  due  to  thaw  consolidation  of  both  the 
mineral  soil  and  the  peat.  There  is  no  evidence  for  cyclical  aggrad¬ 
ation  and  degradation. 

6.3  Open  System  Moisture  Migration  in  Permafrost 

Various  lines  of  evidence  presented  in  this  thesis  indicate  that 
long  term  moisture  migration  through  the  frozen  paisa  soils  may  have 
occurred.  This  flux  has  taken  place  in  spite  of  a  relatively  stable 
frost  front  over  the  past  few  hundred  years.  First,  tritium  values 
in  the  frozen  soils  are  similar  to  those  in  the  sub-permafrost  ground 
water;  their  anomalously  high  values  indicate  mixing  of  modern  water 
with  the  older  permafrost.  Oxygen  isotope  values  in  the  permafrost 
are  best  accounted  for  by  a  conceptual  model  of  regelation  melting 
with  upward  moisture  transfer  into  the  frozen  soils.  Segregated  ice 


* 
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lenses  are  concentrated  in  the  upper  part  of  the  frozen  mineral  soil 
with  massive  zones  of  reticulate  vein  ice  in  the  lower  core;  the  op¬ 
posite  of  what  is  normally  expected  in  aggrading  permafrost  (Mackay, 
1974),  but  consistent  with  a  model  of  post-freezing,  upward  moisture 
flux.  A  single  large  ice  lens  is  usually  found  immediately  above 
the  thawed  zone  in  the  ice  cores,  suggesting  a  continuing  moisture 
supply.  Finally,  the  morphology  of  the  paisa  complex,  with  definite 
marginal  ridges  elevated  up  to  a  meter  above  the  rest  of  the  paisa, 
attests  to  the  sustaining  water  flux  from  the  adjacent  ponded  fen. 

Recently,  laboratory  tests  have  proven  that  water  migration 
occurs  in  frozen  soils  due  to  temperature  induced  potential  gradients 
(Hoekstra,  1966;  Dirksen  and  Miller,  1966;  Ershov  et  al,  1976;  and 
Mageau,  1978).  In  a  series  of  theoretical  papers  Miller  (1972,  1978) 
has  advanced  a  strong  conceptual  argument  that  ice  regelation  is  the 
limiting  process  in  determining  rates  of  water  flux  through  frozen 
soils  containing  segregated  ice. 

The  paisa  complex  environment  is  precisely  the  conditions  under 
which  thermo -migration  of  water  will  most  likely  occur.  The  soils  are 
a  silty  clay  with  potential  unfrozen  water  contents  of  approximately 
30 %;  the  ground  temperatures  are  around  -  1°  G.,  a  nearly  optimum 
temperature  for  unfrozen  water  mobility  (Hoekstra,  1966) .  The  perma¬ 
frost  is  sufficiently  shallow  that  it  is  effected  by  seasonal  temp¬ 
erature  changes  and  there  is  abundant  water  surrounding  the  permafrost 


body. 
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A  semi -quantitative  analysis  can  "be  made  to  determine  whether 
rates  of  thermally  induced  water  movement  are  significant  over  geo¬ 
logic  time.  Ice  movement  hy  the  regelation  process  constitutes  a  vir¬ 
tual  transport  of  sensible  heat  in  the  direction  opposite  to  the  dir¬ 
ection  of  ice  movement  (Miller,  Loch  and  Bresler,  1975) •  The  govern¬ 
ing  equilibrium  mass  and  heat  transfer  equation  can  be  written  (Miller, 
1978): 

-  Ku  dJT  =  -  K  djr  -  ei  L  V  (1) 

d  z  d  z 

where  =  unfrozen  soil  conductivity 

K.p  =  frozen  soil  conductivity 

e^  =  ice  density 


— —  =  thermal  gradient  in  frozen  soil 

L  =  latent  heat  of  fusion 
v.  =  volumetric  ice  flux 

l 

solving  for  the  volumetric  ice  flux: 


vi  u  d  z _ f  d  z  (2) 

e.  L 
1 

parameters  from  the  paisa  soils  are: 

K  =1.2  Cal /  m-hr  °C  K„  =  1.7  Cal/m-hr  °C 

u  * 

=  .05  °c/m  e.  =  910  kgm/  v? 

d  z  1 


Substituting  these  values  into  (2)  yields  a  volumetric  ice  flux  of 

Q  O  2 

approximately  1  x  10  cm/sec  or  0.3  cur/  year-cm  ,  a  rate  comparable 


to  that  experimentally  measured  hy  Mageau  (1978).  Within  the  1500 


year  age  of  the  paisa,  this  represents  a  total 


As  there  is  an  average  thickness  of  1.5  meters  of  segregated  ice  in 
the  paisa,  steady  state  processes  could  theoretically  result  in  a 
redistribution  of  3 0%  of  the  mineral  soil  ground  ice. 

Many  permafrost  bodies  date  from  at  least  the  last  glaciation. 
Ground  ice  in  the  arctic  regions  of  Canada  that  escaped  glaciation 
is  considerably  older.  If  the  steady  state  regelation  rates  calcu¬ 
lated  above  are  representative  values,  then  much  of  the  natural 
ground  ice  must  be  viewed  as  a  dynamic  body  subject  to  continuing 


mass  transfer. 


.  ' 


6  A  Further  Research 


In  North  America,  the  study  of  permafrost  in  disco ntinuously  froz¬ 
en  peatlands  has  not  received  the  same  attention  as  the  more  spectac¬ 
ular  arctic  periglacial  features.  Many  topics  remain  which  require 
further  research.  The  study  described  in  this  thesis  involved  sev¬ 
eral  different  field  and  laboratory  programs,  most  of  which  need  fur¬ 
ther  refinement  or  verification. 

Refinement  in  the  use  of  isotopes  from  peat  as  paleoclimatic  in¬ 
dicators  is  perhaps  the  most  potentially  rewarding  avenue  of  research. 
Acquisition  of  site-specific  temperature  records  would  be  invaluable 
in  deciphering  the  history  of  permafrost  development.  The  ubiquity  of 
peat  in  the  sub-arctic  would  allow  regional  correlations.  If  detailed 
records  were  obtainable,  the  potential  for  mathematically  modeling 
paisa  development  exists.  Those  topics  needing  most  urgent  work,  are: 

1.  an  understanding  of  the  plant  physiology.  Changes  in  the  isotopic 
content  of  Sphagnum  needs  to  be  determined  after  growth  in  a  con¬ 
trolled  environment. 

2.  Further  correlation  of  modern  Sphagnum  isotopes  with  respective 
mean  and  seasonal  growth  temperatures. 

3.  Further  investigation  into  the  role  of  microbial  decomposition  on 
oxygen  isotopes;  possibly  through  the  use  of  nitrogen  isotopes. 

There  seems  little  further  need  for  immediate  work  on  regional 
peatland  development.  Consistent  patterns  seem  to  exist  throughout 
temperate  and  sub-arctic  North  America.  More  productive  research  would 


be  to  analyze  and  date  a  series  of  closely  spaced  cores  across  a  fen- 
pal  sa  sequence,  trying  to  determine  local  plant  successions  and  the 
role  of  vegetative  and  edaphic  factors  in  permafrost  development. 

Both  the  theory  and  testing  of  thaw  consolidation  techniques  are 
well  established  and  require  little  further  laboratory  work.  The  tech¬ 
niques  should  be  used  to  further  test  the  hypothesis  that  peat  consol¬ 
idation  plays  a  major  role  in  paisa  development.  Field  tests  would  be 
particularly  useful,  though  would  require  considerable  imagination  to 
implement. 

The  upward  migration  of  water  through  frozen  soils  is  an  intrique- 
ing  problem.  While  of  academic  interest  for  its  role  in  redistribution 
of  ice  in  permafrost,  the  subject  also  has  considerable  practical  im¬ 
portance.  The  long  term  stability  of  cold-temperature  pipelines  could 
be  affected  by  post-freezing  frost  heave.  While  the  phenomenon  has’ 
been  monitored  in  the  laboratory,  paisa  investigations  could  verify  its 
long  term  occurrence  and  behavior.  Isotopic  techniques  may  be  partic¬ 
ularly  useful  in  this  regard. 
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APPENDIX  Is  Thermistor  Construction 

The  thermistors  were  calibrated,  in  the  lab  at  5  °C.  intervals 
from  -20  °C.  to  +20  °C.  against  a  previously  calibrated  thermistor. 

The  thermistor  leads  were  soldered  at  low  temperatures  to  copper 
wires,  sheathed  in  1.5  mm  'Heat  Shrink'  rubber  tubing,  and  then  in¬ 
serted  into  3  mti  copper  tubing  and  packed  in  liquid  fiberglass.  The 
thermistor  head  was  allowed  to  protrude  .5  mm  out  of  the  copper  tube. 
The  tube  was  sheathed  in  3  mm  Heat  Shrink  tubing,  and  the  copper 
wires  soldered  to  PVC  triple  insulated  dipole  cable.  Several  layers 
of  self -vulcanizing  electrical  tape  were  wrapped  around  the  cable. 
Finally,  the  entire  thermistor  plus  several  centimeters  of  the  cable 
was  dipped  in  'Scothchkote'  several  times. 

The  needle  probe  thermistor  is  initially  prepared  in  a  similar 
way  to  the  permanent  thermistors.  The  thermistor  head  is  then  in¬ 
serted  into  a  pointed,  threaded  brass  rod,  12  cm  in  length.  A  rubber 
insulated  handle  is  screwed  into  the  rod,  the  thermistor  cable 
exiting  at  the  other  end  of  the  handle.  This  device  completely  pro¬ 
tects  the  thermistor  head  from  moisture  and  shock.  The  pointed  rod 
can  be  inserted  into  the  end  of  the  large  diameter  core  immediately 
after  removal  from  the  drill  hole  or  can  free  fall  to  the  bottom  of 
the  hole.  This  allows  more  continuous  temperature  logging,  but  not 
to  the  accuracy  of  permanent  thermistors.  The  needle  probe  thermistor 
was  also  fitted  to  a  flexible  steel  rod  which  was  inserted  to  various 
depths  in  the  unfrozen  bog  and  fen. 
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Five  permanent  thermistors  were  installed  in  the  frozen  paisa  and 
four  in  the  unfrozen  bog.  The  thermistor  cables  were  securely  taped 
to  the  outside  of  the  PVG  pipe  used  for  the  installation  of  piezometers, 
and  lowered  into  the  drill  holes.  Holes  were  back-filled  and  sealed 
in  the  usual  manner. 


* 


Appendix  2 


:  Results  of  chemical  analysis  of  ground  ice  in 
cores  S78,  E3,  and  Ell. 
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A.  Results  of  thaw  consolidation  tests  on  core 
S100 

B.  Results  of  tests  using  empirical  relationships 
to  predict  total  thaw  settlement  in  organic 


and  mineral  soils 
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APPENDIX  4:  Testing  Methods  for  Thaw  Consolidation 

Methods  of  preparing  samples  for  thaw  consolidation  testing  and 
the  apparatus  and  instrumentation  used  in  the  tests  have  been  described 
by  Roggensack  (1978)  in  considerable  detail,  and  will  be  only  briefly 
described  here. 

The  permafrost  oedometer  or  permode  was  the  basic  piece  of  equip¬ 
ment  (Figure  A5) .  It  is  of  the  split  barrel  design  to  facilitate  the 
testing  of  undisturbed  samples.  The  permode  walls  are  constructed  of 
low  thermal  conductivity  PVC  and  lined  with  Teflon  to  eliminate  side 
friction.  De-aired  sintered  stainless  steel  porous  stones  permitted 
drainage  through  the  top  or  bottom  of  the  sample.  The  lower  porous 

stone  was  connected  to  a  pressure  transducer  with  a  resolution  of  ).15 
2 

KN/ni  and  monitored  with  a  Budd  Strain  gauge.  This  was  used  to  monitor 
pore  water  pressures  during  thaw  and  for  constant  head  permeability 
tests.  A  linearly  variable  displacement  transducer  (LVDT)  measured 
sample  settlements  to  an  accuracy  of  better  than  2  x  10  0  meters. 

The  soil  sample  is  turned  on  a  lathe  to  a  diameter  of  .0 64  m  to 
closely  fit  the  permode.  Samples  are  then  milled  to  an  average  length 
of  .07  m,  inserted  into  a  vacuum-greased  rubber  membrane  and  placed  in 
the  permode.  Rubber  0-rings  at  the  load  cap  and  base  pedestal  prevent¬ 
ed  leakage.  Care  was  taken  to  eliminate  any  trapped  air  between  the 
specimen  and  the  loading  cap  and  base.  The  apparatus  is  de-aired,  re¬ 
moved  from  the  cold  room,  and  then  mounted  in  the  loading  frame. 
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Thawing  was  uncontrolled  at  room  temperature  and  under  undrained 
conditions.  After  complete  thawing  at  a  stress  equal  to  approximately 
one-third  the  overburden  stress,  the  permode,  still  in  the  undrained 
condition,  was  sequentially  loaded  with  weights  up  to  the  overburden 
stress.  Maximum  pore  pressures  with  each  load  increment  were  recorded. 
This  technique  allows  the  determination  of  the  residual  stress  of  a 
soil  (see  Roggensack,  1977;  Nixon  and  Morgenstem,  1973)*  The  weights 
were  then  removed,  and  the  undrained  consolidation  recorded.  The 
drains  at  top  and  bottom  of  the  sample  were  opened  and  the  permode  re¬ 
loaded.  Coefficients  of  consolidation  were  determined  in  the  usual 
manner.  The  soil  was  allowed  to  consolidate  for  24  hours  under  each 
load  increment  and  the  relative  settlement  recorded.  The  coefficients 
Aq  and  aQ  were  determined  using  the  technique  outlined  by  Luscher  and 
Afifi  (1973) •  After  total  overburden  stress  had  been  reached  the  sam¬ 
ple  was  allowed  to  consolidate  for  72  hours  and  the  total  thaw  consol¬ 
idation,  E°,  recorded. 


The  overburden  pressures  were  determined  from  the  relationship 


where 


V  p’ 

Vs* 


=  height  of  paisa  above  water  table 

=  depth  of  frozen  peat  beneath  the  water  table 

=  depth  of  mineral  soil  beneath  the  water  table  to  the  thaw 
front 

=  thawed  unit  weight  of  peat 
=  submerged  unit  weight  of  peat 
=  submerged  unit  weight  of  mineral  soil 


APPENDIX  5:  Stefan’s  formula. 


The  simplest  method  for  calculating  frost  penetration  depth, 
X,  in  a  homogeneous  soil  is  Stefan’s  formula: 


X 


Z  K  dT-t 

L 


where:  K  -  coefficient  of  thermal  cond. 
dT  -  step  temperature 
t  -  time 


L  -  latent  heat  of  fusion 


APPENDIX  6:  Testing  Methods  and  Apparatus  for  Frost  Heave 
Characteristics  of  Paisa  Soil 

The  experimental  apparatus  and  methods  used  to  obtain  the  basic 
frost  heave  data  is  identical  to  that  described  by  Hill  and  Morgen- 
stern  (1977)  and  very  similar  to  that  used  by  Penner  and  Ueda  (1977). 
Briefly,  the  test  cell  is  a  PVC  lined,  heavily  insulated  chamber  cap¬ 
able  of  holding  a  pre-consolidated  soil  sample  0.10  m  in  diameter  and 
up  to  .15  m  in  height.  The  sample  is  sheathed  in  a  thin  rubber  mem¬ 
brane  to  prevent  side  friction.  Four  thermocouples  are  imbedded  in 
the  cell  walls.  Consolidation  and  overburden  pressures  were  provided 
by  a  hanging  weight  loading  apparatus  mounted  above  the  freezing  cell 
Freezing  cold  side  temperatures  were  introduced  at  the  upper  plate  by 
a  circulating  ethylene  glycol  bath.  The  lower  plate  was  kept  at 
temperatures  slightly  above  freezing  by  a  circulating  water  bath.  The 
freezing  cell  was  kept  in  a  constant  temperature  (+2  °C.)  cold  room  to 
reduce  thermal  effects  from  the  sides.  The  base  of  the  sample  had 
free  access  to  de-ionized,  de-aired  water  with  the  water  table  main¬ 
tained  at  the  level  of  the  base  of  the  sample. 

After  consolidation  was  complete,  the  pressure  was  reduced  to  a 
selected  value  and  a  freezing  temperature  imposed  at  the  top  of  the 
specimen.  Heaving  was  "open  system"  throughout  the  test. 

Temperatures  at  the  top,  base  and  sides  of  the  sample,  rates  of 

vertical  heave  and  rates  of  water  intake  by  the  freezing  soil  were 

measured  every  hour  by  a  data  acquisition  system  with  printing  and 
plotting  facilities. 
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APPENDIX  7. :  Frost  Heave  Characteristics  of  the  Paisa  Soils 

The  frost  heave  susceptibility  of  the  soil  is  an  extremely  impor¬ 
tant  characteristic  of  northern  terrain.  This  parameter  is  of  signif¬ 
icance  to  geomorphic  process  studies  in  that  it  defines  the  upper  lim¬ 
its  for  the  development  time  of  permafrost  features  and  allows  an 
assessment  of  the  relative  importance  of  the  various  environmental 
parameters  in  the  development  of  a  cryogenic  landscape.  In  this  Ap¬ 
pendix,  the  frost  heave  charac t eristics  of  the  paisa  complex  mineral 
soil  are  determined.  Results  from  these  tests  are  applied  to  the  prob¬ 
lem  of  the  rate  of  development  of  the  ground  ice  profile  in  Chapter 

The  testing  method  involves  placing  a  temperature  gradient  across 
a  frozen  soil  under  a  specified  overburden  pressure.  After  steady 
state  temperature  conditions  in  the  sample  are  achieved,  the  frost  cell 
is  opened  to  drainage  and  frost  heaving  begins.  The  rate  of  heaving 
of  the  soil  is  monitored  under  different  cold  side  temperatures  and 
various  overburden  pressures.  The  temperature  gradient  in  the  frozen 
portion  of  the  soil  is  approximately  linear,  hence  the  position  of  the 
0  °C.  isotherm  is  determined  by  the  cold  side  temperature  (Figure  A8-l). 
Two  sets  of  tests  were  completed;  typical  cumulative  heaving  rates, 
water  intake  rates  and  temperature  profiles  for  the  tests  are  included 
in  Appendix  6.  The  testing  results  have  been  summarized  in  Table  A8-1. 

During  freezing  of  soils  in  the  field,  the  position  of  the  free 
water  table  beneath  the  frost  front  varies  according  to  the  suction 


. 


-consolidation  Pressure  during  Cold  Side  Temp,  gradient 

ssure  (Kg  /  2)  heaving  (Kg  /  2)  temperature  in  frozen  zone 

cm  L  (  C/cm) 
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COLD  SIDE 
TEMPERATURE  (°C) 


-4  -2  0 


PORE  PRESSURE 
( Pascals  x  I05  ) 


COLD  SIDE 
TEMPERATURE  (°C) 


-4  -2  0 


PORE  PRESSURE 
(  Pascals  x  I05  ) 


-  PORE  PRESSURE  PROFILE 

- TEMPERATURE  PROFILE 

0°C  ISOTHERM 


FIGURE  a8-1:  CHANGES  IN  PORE  PRESSURE  GRADIENT 

WITH  CHANGES  IN  COLD  SIDE 
TEMPERATURE  . 


potential,  the  soil  permeability  and  the  rate  of  advance  of  the  freez¬ 
ing  front  (see  McRoberts  and  Nixon,  1973)*  In  "the  frost  permode  how¬ 
ever,  a  stationary  freezing  front,  hence  quasi  steady  state  condition, 
are  achieved  in  a  few  hours.  Throughout  the  tests  the  water  table  is 
held  constant  at  the  base  of  the  sample,  creating  two  water  pressure 
boundary  values:  P  =  -  suction  potential  at  the  freezing  front  and 
P  =  atmospheric  at  the  porous  stone.  The  cold  side  temperature  deter¬ 
mines  the  steady  state  position  of  the  freezing  front,  hence  with  sam¬ 
ples  of  identical  thickness  the  distance,  dX,  over  which  the  hydraulic 
pressure  difference,  dtt,  is  dissipated  is  a  function  of  the  cold  side 
temperature.  Figure A8l  illustrates  the  situation  for  a  cold  side 
temperature  of  -4  °G.  and  of  -1  °G.  from  test  BB-1,  It  is  clear  that 
the  thickness  of  soil  through  which  the  water  must  flow  to  reach  the 
freezing  front  is  considerably  greater  (3.3  x)  for  the  1  °  C.  cold 
side  temperature. 

Suction  pressures  generated  at  the  freezing  front  can  be  calcu¬ 
lated  from  the  data  in  Table  A8-1  using  Darcy's  law: 

dH  =  Q  dX 
kA 

where  k  =  soil  permeability  (1.2  x  10  cm/sec) 

A  =  cross  sectional  area  of  soil  sample 
dX  =  distance  between  the  frost  front  and  the  porous  stone 
The  suction  pressures  are  listed  in  Table  A8-1.  They  are  plotted  as 
a  function  of  cold  side  temperature  in  Figure  A8-2  and  as  a  function 
of  the  temperature  gradient  in  the  frozen  zone  in  Figure  A83 .  The 
plots  are  very  similar  in  form.  There  appears  to  be  a  limiting 
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CALCULATED  SUCTION  AT  THE  FREEZING 
FRONT  AS  A  FUNCTION  OF  COLD  SIDE 
TEMPERATURE  . 


FIGURE  A8-2 


SUCTION  AT  FREEZING  FRONT 
(  Pascals  x  I05  ) 
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TEMPERATURE  GRADIENT  IN 
FROZEN  SOIL  (°C/cm.) 


CALCULATED  SUCTION  AT  FREEZING  FRONT 
AS  A  FUNCTION  OF  TEMPERATURE 
GRADIENT  IN  THE  FROZEN  SOIL  . 


FIGURE  A8-3 


suction  potential  at  some  temperature  between  0  °G,  and  -1  °G.  At 
cold  side  temperatures  below  -1  °C . ,  suction  potentials  are  nearly 
constant.  This  conclusion  agrees  with  those  of  Hill  (1976)  who 
also  found  a  limiting  suction  in  freezing  soils.  At  cold  side 
temperatures  below  -1  °C . ,  pressure  is  the  dominant  variable  in  de¬ 
termining  heave  rates.  The  relationship  derived  by  Linnel  and  Kepler 
(1959)  between  In  heave  rate  and  overburden  pressure  is  supported 
(see  Appendix  9).  In  the  temperature  range  0  0  to  .5  °C.,  the  temp¬ 
erature  variable  is  by  far  the  most  important  variable.  Small  temp¬ 
erature  depressions  at  the  freezing  front  or  small  changes  in  the 
temperature  gradient  of  the  frost  fringe  (Miller,  1978)  produce  large 
increases  in  the  interfacial  suction  potential  in  the  vicinity  of  0  °C. 

Penner  and  Ueda  (1977)  and  Penner  and  Walton  (1978)  have  main¬ 
tained  that  the  hydraulic  potential  (dH)  is  a  direct  function  of  the 
cold  side  temperature.  Results  presented  above  show  that  it  is  the 
hydraulic  gradient  (dH/cLX)  and  not  hydraulic  potential  that  is  a 
function  of  cold  side  temperature.  Further  discussion  of  the  ram¬ 
ifications  of  this  error  in  the  frost  heave  theory  developed  by 
Penner  et  al  is  beyond  the  scope  of  this  thesis. 

In  the  paisa  soils,  a  maximum  overburden  pressure  of  1  kg/  cm 
will  occur.  Approximate  suction  potentials  at  the  freezing  front  of 
3  x  10-5  pascais  occur  at  this  pressure;  this  result  was  used  in  Chap¬ 
ter  5  to  approximate  the  development  time  for  the  observed  ground  ice 
in  the  paisa  complex. 


